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Abstract 

Ultracold neutral plasmas, formed by photoionizing laser-cooled atoms near the 
ionization threshold, have electron temperatures in the 1-1000 kelvin range and ion 
temperatures from tens of millikelvin to a few kelvin. They represent a new frontier 
in the study of neutral plasmas, which traditionally deals with much hotter systems, 
but they also blur the boundaries of plasma, atomic, condensed matter, and low 
temperature physics. Modelling these plasmas challenges computational techniques 
and theories of non-equilibrium systems, so the field has attracted great interest from 
the theoretical and computational physics communities. By varying laser intensities 
and wavelengths it is possible to accurately set the initial plasma density and energy, 
and charged-particle-detection and optical diagnostics allow precise measurements 
for comparison with theoretical predictions. 

Recent experiments using optical probes demonstrated that ions in the plasma 
equilibrate in a strongly coupled fluid phase. Strongly coupled plasmas, in which the 
electrical interaction energy between charged particles exceeds the average kinetic 
energy, reverse the traditional energy hierarchy underlying basic plasma concepts 
such as Debye screening and hydrodynamics. Equilibration in this regime is of par- 
ticular interest because it involves the establishment of spatial correlations between 
particles, and it connects to the physics of the interiors of gas-giant planets and 
inertial confinement fusion devices. 



Contents 

1 Introduction d 

2 Experimental Methods 
2.1 Creation of ultracold neutral plasmas @ 

Preprint submitted to Elsevier 2 February 2008 



2.2 Detection of ultracold neutral plasmas [l(j 

3 Theoretical Description 

3.1 Macroscopic Approaches LL7 

3.2 Microscopic Approaches [27 

4 Physical Processes in Ultracold Neutral Plasmas [33| 

4.1 Initial Electron Equilibration |34 

4.2 Initial Ion Equilibration 

4.3 Collective Electronic Plasma Modes 

4.4 Plasma Expansion 

4.5 Electron Heating Mechanisms 

4.6 Coulomb Coupling Parameters [H3 

5 Achieving Strong Coupling |64 

5.1 Ionic coupling |65 

5.2 Electronic coupling [6^ 

6 Conclusions and Future Directions [71 
A Important Quantities [73| 
B Extracting Plasma Parameters from the Absorption Spectrum [ll 
C The Quasineutral Approximation 
References \7a 



2 



1 Introduction 



In a conventional neutral plasma, ions and electrons are created from atoms 
and molecules by ionizing collisions between particles. Because a typical ion- 
ization potential is on the order of an electronvolt, most neutral plasmas have 
temperatures of thousands of kelvin or more. (See Fig. [TJ) Plasma physics 
in this regime encompasses a wealth of fascinating fundamental phenomena 
such as collective modes, instabilities, and transport mechanisms, and it can 
contribute to important applications such as lighting technologies, plasma pro- 
cessing, and the pursuit of fusion energy. 

Ultracold neutral plasmas [2J stretch the boundaries of traditional neutral 
plasma physics. They are formed by photoionizing laser-cooled atoms near 
the ionization threshold and have electron temperatures ranging from 1-1000 K 
and ion temperatures of around 1 K. While providing access to a new regime, 
they have also proven to be clean and simple systems that form an excellent 
testing ground for basic plasma theory, and they elucidate physics found across 
a much wider spectrum of plasma energy and density. 

The clearest distinguishing feature of ultracold neutral plasmas is that par- 
ticles can be in or near the strongly coupled regime [3J. In strongly coupled 
plasmas the Coulomb interaction energy between charged particles exceeds the 
average kinetic energy, which is quantified by the Coulomb coupling parameter 



where a = [3/(47rp)] 1//3 is the Wigner-Seitz radius characterizing the separa- 
tion between particles in a plasma of density p and temperature T. T > 1 
reverses the traditional energy hierarchy that underlies our normal under- 
standing of plasmas based on concepts such as Debye screening and hydro- 
dynamics. Strongly coupled plasmas exist in dense astrophysical systems [I], 
matter irradiated with intense laser fields [5f6f7f8], colloidal or dusty plas- 
mas of highly charged macroscopic particles [9][T0~] . and non-neutral trapped 
ion plasmas that are laser cooled until they freeze into Wigner crystals |llj . 
Strong coupling is manifested through the presence of spatial correlations be- 
tween particles, which can be important for the equilibration and collective 
dynamics of the system. Optical probes [T2"irT3] demonstrated that ions in an 
ultracold neutral plasma can be strongly coupled. Electrons seem to remain 
just outside of this regime because the rate of recombination of electrons and 
ions to form neutral particles diverges rapidly as the electron temperature 
drops. 

There is great similarity between the dynamics of ultracold neutral plasmas 
and equilibration of plasmas created by fast-pulse laser irradiation of solid 
[TlfTo] and thin film targets [T(3lll7lll8|ll9ll20ll21f22ll23ll24] . rare gas clusters [25] , 
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Fig. 1. Overview of neutral plasmas in the density-temperature parameter plane. 
Lines of constant Coulomb coupling parameter T (see Eq. ([T])) are indicated. After 

&■ 

and gas jets [26|27j . However, the fast-pulse systems have densities close to 
or higher than the condensed phase. This complicates observation of plasma 
dynamics, which takes place on the scale of the inverse plasma frequency, 
Up 1 = [me /(e 2 p)] 1 / 2 , and is of the order of attoseconds for electron and fem- 
toseconds for ions. Ultracold plasmas, on the other hand, reach the strong 
coupling regime despite the fact that they are extremely dilute (Fig. [[]). The 
low density implies that ultracold plasmas evolve on a time scale that is more 
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easily accessible experimentally. Studies of ion equilibration have revealed os- 
cillations of the ion kinetic energy [28|29ll30ll3T] that are a fundamental char- 
acteristic of equilibration of strongly coupled systems. Expansion of a plasma 
into the surrounding vacuum may be directly observed with excellent spatial 
and temporal resolution in ultracold systems [32 33 34, 13J . This phenomenon is 
important in fast-pulse experiments, especially in the context of the generation 
of high energy particles [T4lfl1)lfll^^ inertial- 
confinement fusion experiments [35], and x-ray lasers [36J. 

High energy-density and strongly coupled plasmas are often difficult to de- 
scribe theoretically because of the importance of many-body interactions, and 
these problems are often best addressed with computational techniques. Many 
body interactions are also important in ultracold neutral plasmas, and compu- 
tational physicists have made significant contributions to the understanding 
of these system by applying techniques usually used for high energy-density 
systems [371138] . Ultracold plasmas have motivated the development and ap- 
plication of specially adapted numerical techniques to describe the presence 
of spatial correlations between the ions [39] and the evolution of electrons 
and ions on very different timescales [27J. Using these techniques it has even 
been possible to lead the experimental efforts by predicting the possibility of 
Wigner crystallization in an ultracold expanding neutral plasma |34j and the 
conditions that lead to it. Again, the excellent experimental control and di- 
agnostics available imply that these computational results can be tested with 
great precision in ultracold systems. 

Immediately after the creation of an ultracold neutral plasma through pho- 
toionization, the system is far from equilibrium, and most research in the 
field has focused on the equilibration process. Theoretical and computational 
work [33|37ll38ll4U|l41ll42ll43j has identified processes that heat the plasma. For 
plasmas in or near the strongly coupled regime, where the potential energy is 
comparable to or exceeds the thermal energy, excess potential energy arising 
from a random distribution of ions [40] or electrons [41] is converted to ther- 
mal energy. This is called disorder-induced heating. Another important heating 
mechanism is three-body recombination (TBR), in which an ion and two elec- 
trons collide to form an energetic electron and an excited neutral atom. The 
TBR rate varies with temperature as T~ 9//2 [44], and at low temperatures it 
can become the fastest process in the plasma. TBR is a crucial process for 
the formation of cold antihydrogen through positron-antiproton recombina- 
tion [451T46] . Conditions in these anti-matter plasmas are similar to ultracold 
neutral plasmas except a large magnetic field is usually present to contain the 
particles. Surprisingly, recombination in ultracold neutral plasmas also has 
similarities to stellar dynamics in globular clusters [47] . 

The study of collective modes has proven to be a rich area of research, orig- 
inally motivated by the possibility of observing modification of the mode 
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Fig. 2. Schematic setup of a typical ultracold plasma experiment, exemplified with 
strontium atoms. Neutral atoms are laser cooled and trapped in a magneto-optical 
trap operating on the 1 Sq — 1 Pi transition at 461 nm, as described in [51| . In a 
second step, l P\ atoms are ionized by photons from a laser at ~ 412 nm (see Fig. 
0^). Finally, the ionic plasma component is imaged using the 2 5 , 1 / 2 — 2 Pi/2 transition 
in Sr + at 422 nm (see Fig. [3|3) . Reused with permission from [12J. Copyright 2004, 
American Physical Society. 

structure due to strong coupling [48] . Electron plasma oscillations [32.49J and 
Tonks-Dattner modes [50] have served as probes of plasma dynamics and rep- 
resent interesting fundamental problems in their own right. 

The review is structured such that Sections [2] and [3] provide the experimental 
and theoretical basis that allows one to understand and manipulate ultracold 
neutral plasmas. Section H] of the review describes the physical processes which 
take place in an ultracold neutral plasma from the time it is created until it 
expands into the surrounding vacuum. In Section [5j we discuss proposals and 
prospects for achieving very strong coupling of the plasma - in some sense 
the "holy grail" of ultracold plasma physics. We conclude with Section [6] by 
discussing some possible future directions of this relatively young research 
field. 



2 Experimental Methods 



2. 1 Creation of ultracold neutral plasmas 



The production of an ultracold neutral plasma starts with laser-cooled and 
trapped neutral atoms. There are many good resources on this topic, such as 
the book "Laser Cooling and Trapping," by Metcalf and van der Straten |52j . 
In a configuration of laser beams and magnetic fields known as a magneto- 
optical trap (Fig. [2]), as many as 10 9 atoms can be cooled to millikelvin or even 
microkelvin temperatures, and confined at densities approaching 10 11 cm -3 . 
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Fig. 3. Strontium atomic and ionic energy levels with decay rates relevant for ul- 
tracold plasma experiments. (A) Creation of an ultracold plasma starts with laser 
cooling and trapping neutral atoms using the Sq — l P\ transition at 461 nm. 1 P\ 
atoms are then photoionized with photons from a laser at ~ 412 nm. (B) Ions are 
optically imaged using the 2 Si/2 — 2 P\/2 transition in Sr + at 422 nm. 2 P\/2 ions 
decay to the 2 D 3 / 2 state 7% of the time, after which they cease to interact with the 
imaging beam. This does not affect most experiments because ions typically scatter 
fewer than one photon while the imaging beam is on. The level scheme is similar 
for calcium |13j . 

The typical spatial distribution of the cloud is of spherical Gaussian shape, 

p(r) = p exp(-r 2 /2a 2 ) , (2) 

with a width of o = 200 — 1500 pm. The first ultracold plasma experiments 
were performed at the National Institute of Standards and Technology in 
Gaithersburg (NIST) with metastable xenon [53], but any atom that can be 
laser cooled and photoionized may be used, such as strontium [5Tf5"4"] . calcium 
[13], rubidium [551156] . and cesium [56J. A plasma has also been produced by 
photoionizing a Bose-Einstein condensate [57J. 

A pulsed dye laser excites the laser-cooled atoms above the ionization thresh- 
old (Fig. [3]A.). Because of their light mass, the electrons have an initial kinetic 
energy (E e ) approximately equal to the difference between the photon energy 
and the ionization potential. E c /kB can be as low as the bandwidth of the 
ionizing laser, which is ~ 100 mK with standard pulsed dye lasers, but most 
studies so far have dealt with E e /kB between 1 and 1000 K. The initial kinetic 
energy for the ions is in the millikelvin range. The number of atoms ionized 
(Ni), and thus the density of the plasma (p), is controlled by varying the 
energy of the photoionizing laser pulse. 

We will discuss the dynamics of the plasma in greater detail in Section HJ 
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Fig. 4. Stages and timescales in the dynamics of an ultracold neutral plasma. Elec- 
tron equilibration is a fast process in ultracold neutral plasmas, and it occurs in a 
time equal to the inverse electron plasma frequency, u>~g = yj m e eo/(p e e 2 ). Anal- 
ogously, the timescale for ion equilibration in the second stage is the inverse ion 
plasma frequency, uj~\ = \J m\_eQ/{p\e 2 ). The timescale for plasma expansion is set 
by the hydrodynamic time, r exp = \/micr 2 / ksT e . 

but Fig. H] illustrates how the evolution can be divided into three phases. The 
fastest process is the equilibration of the electrons because they are relatively 
energetic and light. The more sluggish ions then equilibrate with themselves 
on a slightly longer timescale. And finally, the plasma - which is unconfined 
in typical experiments - expands into the surrounding vacuum. 



2.1.1 Photoionization 

The photoionization process is crucial in these experiments, and the book 
"Rydberg Atoms" [58] by T. Gallagher provides an excellent introduction 
to this subject. Most ultracold plasma experiments have used two-photon 
ionization, in which one photon from the continuous-wave (CW) cooling laser 
excites atoms to the upper level of the cooling transition, and a pulsed laser 
then excites them to the continuum. (See Fig. [3J) 

In equilibrium the fraction of atoms in the upper state of the cooling transition 
can approach 50%, and is given by 



= s /2 
J l + So + (25/7) 2 ' {) 

where 5 is the detuning from resonance of the CW laser, 7 is the full linewidth 
at half maximum for the transition, and So = I / h&t is the saturation parameter 
for the excitation beam of intensity I. J sat is typically a few mW/cm 2 , and 
it is the intensity required to drive the transition at a rate equal to 7. The 
excited-state fraction during normal MOT operation is typically about 10%, 
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but an extra excitation beam may also be used to attain / fa 0.5. Highly 
saturating the transition with a pulsed- dye- amplified beam yields even better 
results 1 131. 



Photoionization from the intermediate state is best described in terms of a 
cross section for absorption of ionizing photons, opi. Cross sections near the 
ionization threshold are of the order of the atomic unit of area, opi ~ cZq ~ 
10~ 18 cm 2 . The probability of ionization is given by 1 — exp(— Fapi), where 
the number of photons per unit area in the pulse is F = J pulse dtl / hv^^ for 
instantaneous laser intensity I. A blue (450 nm), 1 mJ, 5 ns pulse has about 
2 x 10 15 photons, and an intensity of 20 MW/cm 2 for a pulse area of 1 mm 2 . 
This gives a 2 x 10 17 cm -2 integrated photon flux and an ionization probability 
of 20% for atoms in the intermediate level. With more power it is possible to 
ionize nearly 100% of the atoms in the intermediate state. 

Cross sections are normally constant over the near-threshold-ionization region 
of interest for ultracold plasmas, unless there is some physics feature that mod- 
ifies them, such as an autoionizing resonance, which is the case in strontium 
|59j . In alkali metals except lithium, cross sections for single-photon excita- 
tion from the S ground state to near the ionization threshold show Cooper 
minima [58J. This is one reason why most ultracold plasma experiments use 
two-photon excitation. 

The ionization process adds momentum and energy to the electrons and ions. 
The particles can be treated as at rest initially, and recoil during excitation to 
the intermediate state is negligible. Momentum and energy conservation for 
the pulsed laser photons of frequency v gives 

2 2 

AE = hu - E IP = A. + , (4) 
2m ; 2m e 

Kk = pi + p e , (5) 

where the ionization potential is Ejp. The photon momentum k = |k| is small. 
For instance, at 10 K, hk 2 ~ (p 2 )/100, which implies that p\ fa p 2 . The ion 
kinetic energy is thus much smaller than the electron kinetic energy, E e . Hence, 
the electrons take away essentially all the excess photon energy, AE, as kinetic 
energy: 

2 

E e = ^fa AE. (6) 



2m e 



The ion kinetic energy is then 



Pi _ Pe 



AE— , (7) 



2mi 2m ; m 

which is only on the order of millikelvin even for AE = Zc^lOOO K. Within a few 
nanoseconds after photoionization [60J, the electrons collisionally thermalize 
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with themselves at ksT c « \E e . We will describe below how processes in the 
plasma can subsequently heat the electrons and ions to substantially higher 
temperatures. 

2.1.2 Spontaneous ionization of a dense Rydberg gas 

An ultracold neutral plasma also forms from spontaneous ionization in a dense 
cloud of ultracold Rydberg atoms. In fact, the first experiments focused on this 
phenomenon [61] , and the hope was to observe a gas-phase analog of the Mott 
insulator-to-conductor transition familiar from condensed matter physics [62] . 
The idea was first explored by Haroche and co-workers [63J when they saw 
spontaneous ionization in a beam of Rydberg atoms. 

In experiments performed so far, collisions, as opposed to a phase transition, 
drive the spontaneous ionization process in beams or laser-cooled clouds, and 
after the formation of the plasma, the dynamics is very similar to the evolu- 
tion of a system created by direct photoionization [33]. We will touch upon 
this work briefly in Sec. 14.1.21 However, it will not be the main focus of this 
review, and for a deeper discussion we refer to papers exploring the sponta- 
neous evolution of Rydberg atoms into a plasma [61] and dynamics of Rydberg 
atoms embedded in an ultracold plasma [55]. Related questions regarding the 
stability of a Rydberg gas are important for fundamental interest and also 
because this is a proposed architecture for quantum computing [65]. 

2.2 Detection of ultracold neutral plasmas 

2.2.1 Charged particle detection 

The most common diagnostics for studying ultracold plasmas is detection of 
charged particles after they have left the plasma. The basic idea is shown in 
Fig. [5J The initial experiments at NIST on plasmas produced by direct pho- 
toionization [2,32,66,67] detected electrons. Experiments exploring the spon- 
taneous evolution of Rydberg atoms into a plasma [61] and dynamics of Ryd- 
berg atoms embedded in an ultracold plasma [55|68] have used both ion and 
electron detection. 

For monitoring electrons that escape the plasma during normal evolution, it 
is essential to not perturb the plasma with large electric fields. However, one 
needs a field of about 5 mV/cm to guide the escaping charged particles to the 
detector. This produces a time-of-flight delay for electrons until they reach 
the detector of about 1 /xs. Since electrons move away from the plasma in 
all directions, it also means that the detection efficiency is not 100%. easily 
be ramped to much larger values that can quickly dump all charged particles 
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Fig. 5. Charged particle detection diagnostics. When electrons escape from the 
Coulomb well formed by the ions, electric fields created by potentials on wire-mesh 
grids direct them to an electron multiplier or microchannel plate. By reversing the 
electric field, ions can be detected. The grids can also be used to apply oscillating 
radio frequency fields to excite plasma collective modes. 

of a given sign on the detector. With careful design, values well in excess of 
100 V/cm can be obtained, which are also useful for pulsed field ionization of 
Rydberg atoms [58J. 



2.2.2 Optical probes 

Probing an ultracold plasma through charged particles has some limitations. 
For example, remote detection of particles tends to only provide information 
on average properties of the plasma, and cannot resolve variations, such as 
ion acoustic waves or ion-ion spatial correlations. The time resolution is also 
limited by the time of flight for electrons to reach the detector. 

For detecting electrons, there are currently no options besides charged parti- 
cle diagnostics. For probing ions, however, optical methods offer a powerful 
alternative that can provide in situ, non-destructive measurements, with ex- 
cellent spatial, temporal, and spectral resolution. Absorption imaging [12] has 
been used to probe ion-ion equilibration and expansion of the plasma during 
the first few microseconds after photoionization. Fluorescence monitoring has 
given information on the plasma expansion dynamics [13] from shortly after 
photoionization to 50 /xs later. Both techniques have great potential to study 
phenomena such as ion collective modes [69], shock waves [33J, recombination, 
and particle-particle spatial correlations |llj . 

Absorption imaging is an adaptation of one of the most powerful techniques 
for studying laser cooled and trapped neutral atoms [52]. Implementation of 
this probe to study ultracold plasmas requires much higher temporal resolu- 
tion than typically needed for neutral atom experiments because the plasma 
evolves quickly - in as short as nanoseconds in some cases. It is also essential 
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Fig. 6. Optical depth of a strontium ultracold neutral plasma. The delay between 
the formation of the plasma and image exposure is 85 ns. The plasma contains 
7 x 10 7 ions and the initial peak ion density is poi = 2 x 10 10 cm -3 . Reused with 
permission from [12J. Copyright 2004, American Physical Society. 

to perform experiments with an ion whose ground state possesses electric- 
dipole allowed transitions in the optical regime where cameras and lasers are 
available. Alkaline-earth elements possess this property, which is why they are 
used in quantum computing and optical frequency standard experiments us- 
ing trapped ions. Strontium is a good choice for plasma experiments because 
its laser cooling is well-developed [51 J and imaging can be performed using 
the Sr + 2 S\/2 — > 2 P\/2 transition at 422 nm (Fig. EB)- Calcium is also a good 
choice because it has a similar level structure [T3|70] . 

To record an absorption image of the plasma, a collimated laser beam, tuned 
near resonance with the principal transition in the ions (Fig. EB), illuminates 
the plasma and falls on an image-intensified CCD camera. Ions scatter photons 
out of the laser beam and create a shadow that is recorded by the camera. 
The optical depth (OD) is defined in terms of the image intensity without the 
plasma (J Q ) and with the plasma present (I) as 



OD(x,y) = ln[I (x,y)/I(x,y)} . (8) 

Figure M shows a typical absorption image. By varying the delay At between 
the formation of the plasma and image exposure, the time-evolution of the 
plasma can be studied. The minimum camera exposure gate width for standard 
intensified CCD cameras is ~ 10 ns. 

Valuable information is provided by the absorption spectrum, which can be 
recorded by varying the frequency of the probe beam v [12]. For an average 
spectrum of the entire plasma, the optical depth is integrated over the plane 
perpendicular to the laser, 
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Fig. 7. Absorption spectra of ultracold neutral plasmas (Eq. Q). The frequency is 
with respect to a Doppler-free absorption feature in a strontium discharge cell. 
Both spectra correspond to T e = 56 K and an initial peak plasma density of 
poi = 2 x 10 10 cm -3 . Data are fit with Voigt profiles (Eq. (|12p ). Note the increase in 
linewidth for longer delay At (see Eq. (|B.4j) in appendix[B]). Reused with permission 
from [54J. Copyright 2005, Institute of Physics. 



In reality the camera pixels are summed, and typical averaged absorption 
spectra are shown in Fig. [7J 

In order to quantitatively interpret the spectra, it is important to account 
for the inhomogeneous density distribution of the plasma and other effects 
that will be discussed in subsequent sections. Among them are thermal ion 
motion, plasma expansion, and the lack of global thermal equilibrium for the 
ions. Detailed descriptions of the spectral analysis were given in [51] and |71j . 
Further aspects are also discussed in App. [B] The OD can be related to 
underlying physical parameters according to Beer's law, 



where pi(x,y,z) is the ion density, and a[u, T^(r)] is the ion absorption cross 
section at the image beam frequency, v. The absorption cross section is a 
function of ion kinetic energy, or temperature, due to Doppler broadening. 
The ion temperature varies with density, so for a spherically symmetric but 
inhomogenous plasma (e.g., for a gaussian plasma cloud as in Eq. T\ — 
T[(r) and therefore a varies with position. Equation ([TO]) can be used to relate 
the spectrum to the density weighted average of the absorption cross section 




(9) 




(10) 
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The absorption cross section a[u, 71 (r)] for ions whose motion can be described 
by a temperature 71 is given by the Voigt profile (see appendix [Bl Eq. (IB.lj) ). 
When this expression, generalized to include the effect of the plasma expansion 
(Eq. (IB. 3D ) is substituted into the expression for the spectrum (Eq. (ITTj) ). the 
result is reasonably complicated. However, to a good approximation [54"f7T] . 
the density-averaged spectrum can be fit by a single Voigt profile character- 
ized by an effective ion temperature (7i, e ff) that describes all the Doppler 
broadening due to ion motion (For the derivation, see appendix [B]) , 



b\v)—Ni / as tj ; — exp 

K) 27T J eS J 1 + 4 p^\ 2 V^M?leff) 



cff ) 



\2 



(12) 



where <Jd(T) = \Jk B T /m-J \ is the Doppler width, and 7 e fj = 7o + 7ins is the 
effective Lorentzian linewidth due to the natural linewidth 7 of the transition 
and any instrumental linewidth 7 ins . The center frequency of the transition is 
uq = c/A, and the "three-star" symbol [72] is a numerical factor that accounts 
for the polarization state of the ions and the imaging light. 

With reasonable assumptions valid for a wide range of plasma parameters, 
T i eff can be related to the average ion temperature in the plasma, T i ave = 
/ d 3 rTi(r), but the relation varies as the plasma evolves. During the first 
microsecond after photoionization and while disorder-induced heating (Sec. 
14.2. ip takes place, the plasma is not in local thermal equilibrium (LTE). (See 
Section H~2l ) Strictly speaking, ion temperature is not a well defined concept 
at in this situation. However, 71 i6 g can always be related to the average ion 
kinetic energy and the rms ion velocity along the laser beam through feeT]^ = 

^(-^kin) — m i u z,rms- 

Once LTE is established, it is possible to assign an ion temperature reflecting 
thermal motion, Tj(r), but it varies with position because the plasma lacks 
global thermal equilibrium [73f7T] . The approach to global thermal equilibrium 
is governed by a heat diffusion equation. This leads to a equilibration time 
equal to 

Tglobal ~ CT 2 /D th , (13) 

where the heat diffusion coefficient for a strongly coupled plasma is predicted 
to be on the order of [71] 

Ah ~ a 2 cj p>i . (14) 

This depends on the plasma density through the Wigner-Seitz radius and the 
ion plasma oscillation frequency. r g i b a i is typically greater than 100 /is, which 
is longer than experiments or simulations have probed. 



The average ion temperature just after establishment of LTE, T ; 



LTE 
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important characteristic scale for thermal ion energy. For times much less 
than 



^exp 



rpLTE 

i,ave 

~ 7V 



rpLTE 

i> ave (A £\\ 



the Doppler shift due to plasma expansion is much less than the Doppler width 
due to thermal ion motion, and the expansion can be neglected. Here, r ex p is 
the characteristic time scale of plasma expansion which emerges from the 
kinetic theory of ultracold plasmas. (See Eq. ff25l) in Section l3.1.1l ) References 
[51] and [7T] show through numerical simulation that under these conditions, 
T i>eff « CT i>ave , where C = 0.95 ± 0.05. 

At later times, the expansion is not negligible. But for a broad range of plasma 
initial conditions, the expansion is quite simple and can be described by as- 
suming quasi-neutrality and a self-similar expansion of the plasma (Section 
I3.1.ip . The effect of the expansion velocity on the Doppler width can then be 
incorporated into the effective temperature of Eq. (1121) as shown in appendix 



Regardless of the exact spatial dependence of 71 (r), the integral of the spec- 
trum for the entire plasma is proportional to the number of ions Ni in the 
plasma and the natural linewidth of the transition. This provides an absolute 
calibration of the signal through 



/ duS{y)= I dv d\ pi^a^T^r)] = 3 * A Na ° . (16) 

If Eq. (|9]) is modified so that the OD is integrated over a subsection of the 
plasma image, then T i<e g describes ion motion in the region bounded by the 
integration limits and extending along the propagation direction of the imaging 
laser. The integral over the spectrum is then proportional to the number of 
ions in that region. Analysis of spectra in annular regions was described in 
[28|7T] and proved very useful for studying kinetic energy oscillations (Section 

In fluorescence experiments, which have been performed with calcium [T5|70] . 
fluorescence from a cylindrical region formed by the intersection of the plasma 
with a tightly focused laser beam is collected on a photomultiplier tube. (See 
Fig. [HJ) This approach has the advantage of recording the full time evolution 
of the system with each photoionizing pulse, but repeated scans must be taken 
with the laser at different positions in order to map the spatial distribution. 

When the only Doppler-broadening of the optical transition arises from ther- 
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Fig. 8. Fluorescence detection in an expanding ultracold neutral calcium plasma. 
Left - The plasma is elongated along z. Regions 1-3 occur sequentially in time and 
show how the plasma expands predominantly in the x-y plane. The laser excitation 
probe propagates parallel to the z axis. Right - Electronic levels of the calcium ion 
important for the experiment. Reused with permission from [70| . Copyright 2005, 
American Institute of Physics. 

mal motion of the ions, and the ion velocity distribution along the laser beam 
can be written as a Maxwellian of the form P(v z ) = exp (—v^/2v^ h ) / \/2ivvth, 
the amplitude of the fluorescence signal on resonance is related to the thermal 
velocity, t> t h, through [13] 

where b = 7 A/2 is the velocity that corresponds to a Doppler shift equal to 
the natural linewidth of the transition, 7. If the temperature of the ions is 
known at some time, the signal can be calibrated and the evolution of the 
ion temperature can be extracted p2]. When Doppler broadening due to the 
expansion velocity of the ions becomes important, an analogous expression can 
be derived that allows one to obtain information on the expansion dynamics 

123 ■ 



3 Theoretical Description 



In the last decades, a variety of methods have been developed to describe equi- 
librium states [75ll76ll77f 781179] . relaxation processes [80]l8Tj . transport proper- 
ties {821183] . and atomic processes [E4 : 85 : 86 ; 87] of plasmas. These methods 
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cover a wide parameter space, ranging from cold, low-density plasmas in gas- 
discharges to warm dense matter systems. A typical ultracold plasma behaves 
classically, which facilitates its description. Indeed, the low particle density in 
an ultracold plasma implies a Bruckner parameter r s = (3/47rp e ) 1 / 3 /a 3> 1 
and an electronic Fermi energy Ep <C ksTg. Hence, the values of both of these 
parameters justify a classical treatment of the free charges. Moreover, due to 
the low electron temperatures, recombination forms atoms in highly excited 
states, such that electron-atom collisions jH] and radiative transitions [88] 
can be described classically as well. One might be tempted to conclude that 
a macroscopic treatment based on kinetic theory incorporating some collision 
phenomena may suffice to describe ultracold plasmas. 

However, as it has turned out this is not true. Even a classical description 
of ultracold neutral plasmas is far from trivial and poses several challenges. 
First of all (and this is of course one reason why ultracold plasmas are attrac- 
tive) the ionic plasma component may be strongly coupled, rendering common 
kinetic theories for weakly coupled plasmas inapplicable. Strongly coupled 
problems are usually tackled numerically with a molecular-dynamics (MD) 
approach. Due to the open boundary conditions and long evolution times of 
ultracold plasmas, however, a full MD description of the plasma dynamics 
constitutes an intractable computational task currently. Moreover, the highly 
non-equilibrium plasma dynamics exhibits different relaxation processes on 
very different time scales, which all have to be appropriately taken into ac- 
count for a reliable and realistic understanding of an ultracold plasma. 

Hence, one may ask if it is possible at all to describe ultracold plasmas on 
a macroscopic level based on kinetic and hydrodynamic formulations. This is 
indeed possible to some extent, as we will detail in Section I3~T1 Starting from 
a collisionless description of the plasma, additional effects such as inelastic 
collisions and strong ion correlations are successively incorporated into the 
model. More sophisticated, but also technically more demanding, approaches, 
which describe the plasma dynamics on a microscopic level, are introduced 
in Section 13.21 They can be used and are necessary to assess the validity and 
limitations of the macroscopic descriptions. 



3.1 Macroscopic Approaches 

3.1.1 Collisionless plasma dynamics 

If we assume the Coulomb coupling parameter of both the electrons and the 
ions to be much less than unity and neglect any collisional processes in the 
plasma, the dynamical equations describing the one-particle phase space den- 



17 



sity f a (r Q , v, 



a 



) [89] simplify to the well known Vlasov equations 



Q 



+ v, 



a 




Of 



a 



dr, 



— m 



a 



Here, a = e, i for electrons and ions 



and charge of species a, and <p is the total mean-field potential determined by 
the Poisson equation 



with p a = I fadv a . As a paradigmatic problem of plasma physics, the colli- 
sionless plasma expansion into vacuum has been considered for a long time 
[90f9~T] and solutions of Eqs. (fT8l) have been investigated in a number of recent 
publications [92|93] . Ultracold plasmas for the first time allow for a combined 
theoretical and experimental study of this problem. Clearly, for the present 
situation of potentially strongly coupled plasmas, the Vlasov equations alone 
will not suffice. But they will serve as the basis of our hydrodynamic approach. 
A suitable approximate treatment of correlation effects will subsequently be 
added step by step, to clarify the influence of correlations on the expansion 
dynamics. 

Due to its nonlinearity a general, closed-form analytical solution of Eqs. (1181) 
cannot be given. Here we will employ the additional assumption of quasineu- 
trality [91], under which a large set of self-similar solutions exists. Of particular 
interest in the context of ultracold plasmas is the case of a Gaussian phase 
space density (see Eq. (j2"2l ). discussed in [92], since this is the typical initial 
shape of the plasma created in a magneto-optical trap. For a spherically sym- 
metric plasma it describes a local equilibrium in velocity space, which holds for 
arbitrarily long times of the plasma evolution. We will use a further approxi- 
mation, enforcing a local equilibrium for the electrons (see Eq.f l20p ). in order 
to treat more general cases of strongly coupled or non-symmetric plasmas as 
discussed below. Due to the very small mass ratio m e /mi, the timescale for 
equilibration of the electrons is much smaller than the timescale for motion 
of the ions. Exploiting this fact, one may safely use an adiabatic approxima- 
tion for the electrons, assuming instant equilibration, such that the electronic 
phase space density is given by a quasistationary distribution 



where the time dependence of the velocity distribution function 0( qs ) is only 
implicit in the time dependence of T e . 

1 Here and in the remainder of this section, greek indices denote particle species 
while roman indices label individual particles of a given species. 




(19) 



f e (r,v,t) = fW (r,v,T e (t)) 

= p e (r, t) (v, T e (t)) oc p c (r, t) exp 
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Substituting this expression into the electronic Vlasov equation ([TBI) together 
with the quasineutrality condition (p c (r) « p;(r)) allows one to express the 
mean-field potential in terms of the ion density 

9(fi . - _i<9p e _ 1 9p i 
e^- = KBT e p e — w /cbTcP; — , (21) 
ot or or 



which together with the ionic Vlasov equation ( fl8l) yields a closed set of equa- 
tions for the ionic phase space density Note, that the condition p e (r) ~ Pi(r) 
does not imply the space charge potential to be zero. Rather, it replaces the 
Poisson equation (Tl9l) by Eq. (T2~T]) . In appendix [C] we provide a more detailed 
account of this fact, which may seem somewhat contradictorily and is rarely 
discussed in the literature. 

For a Gaussian spatial plasma density, Eq. (|2~T1) implies a linear force such 
that the simple ansatz 

— (-?4)-(-? S5 ^ £ ) (22) 

provides a selfsimilar solution of the ionic kinetic equation. Here the index 
k = x,y,z labels the (cartesian) components, a is the rms-radius of the Gaus- 
sian spatial distribution, and 7 is a parameter determining the local mean of 
the velocity distribution, i.e. the hydrodynamic ion velocity. Coincidentally, 
the Gaussian form (122]) corresponds to the initial distribution of ions, which 
reflects the Gaussian density profile and zero mean velocity (i.e. 7 = 0) of 
atoms in the magneto-optical trap. Hence, within the collisionless quasineu- 
tral approximation, the expanding plasma retains its shape for all times, and 
its evolution can be parameterized by the macroscopic parameters cr^, jk, 
and T e only. 

Considering a spherically symmetric plasma of Gaussian form as in Eq. (|2l 
and substituting the ansatz Eq. (|2"2"|) into Eq. (TT5j) yields the following set of 
equations for the plasma parameters 



JU 2 = 2 7 a 2 , (23a) 



8_ _ (k B T e + k B T { ) , 
dt m\a 2 



i= VD e 1 7 -r, (23b) 

01 

d 



9£ -(feT i ) = -2 7 fe B T i , (23c) 

|-(A;BT c ) = -27A ; BT e . (23d) 
ot 

In addition to the conserved total energy E tot = |iY i /c B (T , e + T{) + ^N-jn-^a 2 
Eqs. [22] possess two integrals of motion a 2 T e = const, and o 2 T y = const., re- 
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fleeting the adiabatic cooling of both the electrons and the ions during the 
plasma expansion. Using these constants of motion, a simple analytical solu- 
tion of Eqs. (I2"3"j) describing the collisionless plasma expansion can be found 



a 2 (t) = a 2 (0)(l+t 2 /r c 2 xp ) , 
7(*) 



cxp 



Ti(t) 



1 + t 2 /rl p ' 
1 + ' 



cxp 



r.(o) 



(24a) 
(24b) 

(24c) 

(24d) 



where the characteristic plasma expansion time is given by 



T, 



cxp 



mi<j 2 (0) 
\ k B T c (0) 



and the expansion velocity is 



u 



(r,t) = 7 (t)r . 



(25) 



(26) 



Measurements of the average density dynamics for ultracold plasmas are in 
very good agreement with these simple relations for high initial electron tem- 
perature (T e (0) > 70 K) [32]. However, the same experiments have revealed 
significant deviations from this average density dynamics for lower electron 
temperature. As discussed in [42], they mainly arise from inelastic collisions, 
i.e., from the formation of Rydberg atoms during the plasma expansion. This 
is one consequence of collisions in the plasma which we will discuss next. 



3.1.2 Collisional processes 

There are two types of collisional processes in plasmas: elastic collisions, which 
do not affect the total system temperature T e + TJ, and inelastic collisions, 
which may drastically change the temperature of the electrons. On the short- 
est timescale of several ns, elastic electron-electron collisions lead to an equi- 
libration of the electron velocities, giving rise to the electronic temperature 
of the plasma. Implicitly, we have already accounted for this process in the 
previous section within the adiabatic approximation Eq. f[2"Uj) . which assumes 
an instantaneous relaxation towards a Maxwellian electron velocity distribu- 
tion. This assumption is well justified since the timescale for electron-electron 
collisions is by far the shortest one for a typical ultracold plasma. The ion 
velocity distribution, on the other hand, has the initial Maxwellian form of 
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the atomic velocity distribution. Hence, the ion-ion collision integral vanishes 
and does not change the functional velocity dependence of the ion distribution 
as described by Eq. fT2"2l . 

While electron-electron and ion-ion collisions lead to an individual electron 
and ion temperature, respectively, electron-ion collisions tend to equilibrate 
the ion and the electron temperature. However, elastic electron-ion collisions 
may safely be neglected for the description of a freely expanding plasma[il, 
since the corresponding ion heating rate [60] is much less than the rate of the 
initial disorder induced ion heating to be discussed in Section 14.2.11 During 
later times, the expansion leads to an adiabatic cooling of the ions, but also to 
a reduction of the collisional heating rate. Hence, elastic electron-ion collisions 
remain negligible at all times [94J. 

This is not the case, however, for inelastic electron-ion collisions due to the 
low temperature of the electrons, which leads to efficient recombination into 
bound Rydberg states. The formation of bound states during the evolution 
of the plasma can be described most accurately within a classical molecular 
dynamics treatment taking into account all electron-ion interactions [37J3S]. 
However, as discussed above, such an approach is limited to short times due 
to the fact that the small timescale on which the electronic motion takes place 
has to be resolved. Within the present macroscopic approach the problem 
is simplified significantly by adopting a chemical picture, i.e., by introducing 
neutral atoms as an additional particle species. Naturally, this transition is ac- 
complished by splitting the electron-ion two-particle distribution function f ei 
into a free and a bound part [84] , where the bound part is identified with the 
distribution function / a of the atoms. The evolution equation for the atom dis- 
tribution function then involves an integral over the three-particle distribution 
function describing correlations between pairs of electrons and an ion, which 
accounts for inelastic collisions, i.e. recombination, ionization and electron 
impact (de)-excitation. Without going into detail, the corresponding inelas- 
tic processes can be identified as three-body recombination, electron-impact 
ionization, and electron impact excitation or deexcitation. Consequently, the 
electron-ion kinetic equation ffTSl) splits into separate equations for free elec- 
trons, ions, and atoms. Neglecting atom-atom interactions and still retaining 
the Vlasov approximation, the kinetic equations for ions and atoms read 



dfi dfi dfi dip (3) 
1^7 + Vi ^ em ^lT 7T = J ei ' 

Ot OTi (TV; OT; 



9 fa. 9fa, T /n^7\ 

"aT + V a^— = ^ae • (27) 



2 Note that the situation is drastically different in laser-cooled plasmas considered 
in Section 15.1.11 
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The collision terms on the right-hand sides of Eq. ff271) are given by 
J S = K ion{n, Pe(ri), T e )f a (n, p i; r ; ) - ^ K thr (n, p e (ri), T e )/i(pi, r ; ) (28) 

n n 

and 

^ae = $3 [-^bb(p, Pe(r a ), T e )f & (p, p a , r a ) - if bb (n,p, p e (r a ), T e )/ a (n, p a , r a )] 
+i£tbr(™, Pc(r a ), T c )/i(p a , 

I"a) -^ion 

(n, p e (r a ), T e )f a (n, p a , r a ) , (29) 

where the rate coefficients -ft'tbr; -^ion and i^bb for three-body recombination, 
electron impact ionization and electron impact induced bound-bound transi- 
tions depend on the atomic principal quantum numbers n and p, the electronic 
density p e , and the temperature T e . Since these coefficients describe three- 
body processes, no exact analytical expressions for the rates are known. For 
the highly excited Rydberg states in ultracold plasmas an accurate descrip- 
tion is expected from the rate coefficients of [H], which have been obtained by 
fitting analytical functions to the results of classical trajectory Monte Carlo 
calculations. Since the collision integrals Eqs. (1281) and (1291) introduce a non- 
linearity in the spatial particle densities, an exact selfsimilar solution of the 
corresponding kinetic equations can no longer be found. Nevertheless, we as- 
sume the phase space distribution to be the same for atoms and for ions and 
to still be given by Eq. (12"2"1) . which was shown to yield a good description 
of freely expanding ultracold plasmas by comparison with more sophisticated 
calculations [39J. Based on this approximation and from the moments of the 
distribution function 

(V 2 ) = 3cx 2 = iVf 1 J r 2 /i( ri , v b t) dr^V; (30a) 
(vr) = 3 7 a 2 = Nr 1 J v^/ifo, v i; t) dr^vs (30b) 

(^ 2 ) = 3 + 7 V) = Nr 1 J v?Mv h v u t) dridvi . (30c) 

we finally obtain simple evolution equations for the parameters of the plasma 
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Nj k B T c + fc B Ti 
A^i + iV a m i( j 2 



(31a) 







2 



(31b) 



_ 7 



di 
d_ 

di 



(k B Ti) 



P+ri 



[Rbb{p,n)Af a {p) - R hh {n,p)N a (n)} 



(31c) 




tot 



|iVi [k B T e + k B T, + m i7 V] - XXW^ = const. , (31e) 



where 1Z is the Rydberg constant which has the value 1Z = 13.6 eV in electron 
Volts and 1Z = 2.18 x 10~ 18 J in Si-units. The number of atoms with principal 
quantum number n is denoted by J\f a (n), where J2nN a (n) = N a , and the 
density averaged collision rates are defined as 



As can be seen from the energy relation Eq. (13 left , inelastic collisions not only 
change the total number of atoms and ions, i.e., the degree of ionization of the 
plasma, but also affect the temperature of the electronic component. Three- 
body recombination heats the electron plasma, since one electron carries away 
the excess energy gained in the recombination of another electron with an ion. 
Subsequent de-exciting electron- Rydberg atom collisions further heat up the 
plasma, electron impact ionization and excitation tend to cool the electrons. 

The total three-body recombination rate has a strong temperature dependence 
and scales as J2 n Rtbr{ n ) ^ T~ 9 / 2 . Hence, for very high initial electron tem- 
peratures inelastic collisions are of minor importance and Eqs. (1311) reduce 
to Eqs. (1231) . such that the plasma expansion is well described by the simple 
relations Eqs. (1241) . On the other hand, at low initial T e , when the electronic 
Coulomb coupling parameter approaches unity, three-body recombination and 
subsequent electron-atom collisions quickly heat up the electron gas, strongly 
decreasing T e [HS]- This justifies our assumption of an ideal electron gas, in- 
herent in the kinetic description. The ionic temperature, on the other hand, 
remains practically unaffected by inelastic collisions. Hence, ionic correlation 
effects may not necessarily be negligible, and will be accounted for in the next 
section. 




(32) 
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3.1.3 Strongly coupled ions 



Kinetic theories of non-ideal plasmas are commonly based on a small-r expan- 
sion. For ultracold plasmas such approaches are clearly not applicable, since 
the ion component can exhibit Coulomb coupling parameters on the order 
of or even much larger than unity. Alternatively, one can combine a hydro- 
dynamical treatment of the expansion dynamics with numerical results for 
homogeneous strongly coupled plasmas [39]. The starting point is the kinetic 
equation for the ions, which accounts exactly for ion-ion correlations, 

( d d e dip d \ . , , d f ( dVa \ . , , , 

(33) 

where cy is the ion-ion correlation function. To simplify notation we neglect 
inelastic collisions whose contributions can straightforwardly be added, as de- 
scribed in the previous section. 

As before, from Eq. (133]) we derive the time evolution of the moments of j\ up 
to second order, which together with the ansatz Eq. (122]) yields 



J^x 2 = 2 7 a 2 (34a) 

3 d o 3 9 9 ^7 m ^ 7 m 

2 mi di' lcr = 2 mn ° + 2 + 2 



2Ni 



— J p i (r i )r i F ii (r i ) dr, (34b) 



~ (M + m i7 V) = 7 fc B T c - ^ . (34c) 

The correlation pressure force F i; and the correlation energy Ua are obtained 
from 

3 V pi dpi J dr { 



and 



respectively, where 



Ua = — j p i (r i )ii ii (r i ) dri , (36) 



2 J y 

is the mean correlation energy per particle and gn(r, p(r{)) is the spatial pair 
correlation function of a homogeneous plasma of density p(r{). Using 

^ | Pi(r)rF u (r) dr = J u iiPi dr = -U Vl , (38) 
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Eqs. (I3"4"|) are rewritten as 



da 2 9 

— = 2 7 a 2 (39a) 

^c + ^ + foi 2 (39b) 
ot rriicr 2 

— = - Wi -- 7 [/,-- ¥ (39C) 
^ = -2 7 fc B T e , (39d) 
where the last equation follows from energy conservation, 

771; / 2 \ ,3 



E tot = ^iy {v 2 ) + -N-M + NiUz . (40) 

While the description of the influence of ionic correlations on the plasma 
dynamics via the single macroscopic quantity [/„ constitutes a significant sim- 
plification of the problem, the system of eqs. ( 1391) is still not complete since 
it does not contain an equation for the time evolution of An accurate de- 
scription of the evolution of Ua on a kinetic level is rather complicated, and 
beyond the capabilities of the simple approach pursued here. A useful and rea- 
sonable approximation for the evolution of U$ is the so-called correlation-time 
approximation |96j 

BU-. 77.. - 77.( eq) 

CHJ^ ^ _ U a U u (41) 

Ol Tcorr 

where r corr = ^m^o/ (e 2 p ; ) = u~l is the timescale on which pair correlations 
relax towards their equilibrium value, 

U!T>= / Pi (r)4 Cq) (r)rfr, (42) 



and Ujf is the correlation energy of a homogeneous one-component plasma 
(OCP) of density Pi(r) in thermodynamical equilibrium [97J. It is a special 
property of the Coulomb potential that the temperature scaled correlation 
energy is uniquely determined by the Coulomb coupling parameter r i; such 
that numerical values of Ua can be tabulated over a broad range of IYs 



and accurate formulae exist. The following formula interpolates between the 
low-r Abe limit and the high-r behavior [100J 

= fer ' r3/2 (v^f + iff) • < 43 > 

with A x = -0.9052, A 2 = 0.6322 and A 3 = -y/3/2 - A x j 'yfA~ 2 . 

This expression is valid only for an OCP, i.e., a system of ions embedded 
in a neutralizing negative electron background. Therefore, ultracold plasma 
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experiments provide an ideal testing ground for the validity of the widely used 
OCP model in real two-component systems. While the OCP model should 
yield a good description for r e — > 0, we would certainly expect deviations 
from Eq. (143 p due to electron screening of the ion-ion interaction for low 
initial electron temperatures. For correlated but still weakly coupled electron 
components this effect is well described within the Debye-Hiickel theory by an 
exponentially screened Yukawa-type interaction potential 

2 „-r/A D 

V* = ^- , (44) 

where 



A D = ^s k B T e /(e 2 Pc ) (45) 

is the Debye screening length, characterizing the mean distance on which local 
charge imbalances are screened inside the plasma. As this screening effect 
introduces an additional length scale in the plasma the equilibrium correlation 
energy now depends on two parameters, T\ and k = a/ \ D , 

Ut ] = AfcUr, (U + 1) . (46) 



The quantity U (F- x , k) has been obtained from molecular dynamics simulations 
and tabulated over a wide range of Fi and k values |fUllf l02j. 

The range of applicability of the Debye Hiickel is not strictly defined. The 
theory is rigorously valid only for a < Ad. Ultracold plasma experiments, 
however, showed that Eq. pS"]) yields a surprisingly good description even for 
a Ad (r c 1/3) - i.e. with only one electron per Debye sphere [12] . 

Eq. ( 1411 and Eqs. ( 1391) constitute a closed set of equations, describing the time 
evolution of a neutral plasma of electrons and ions. Compared to more involved 
approaches, such as, e.g., Molecular Dynamics methods, the macroscopic evo- 
lution equations are simple enough to allow for significant physical insight 
into effects of ion correlation on the plasma dynamics. While Eq. (139dl) shows 
no direct influence on the electron temperature, the time evolution of the ion 
temperature is significantly modified by the interionic correlations. In addition 
to the adiabatic cooling, which is given by the first term on the right-hand side 
(rhs) of (I39cp in complete analogy to the electronic component, the ionic tem- 
perature is increased (note that Ua < 0) by the development of correlations 
(last term on the rhs of (I39cl) ). This effect, which has been named "disorder- 
induced heating" or "correlation-induced heating," has been widely studied 
for homogeneous systems [80ll37ll41ll40f 103II104II73II30] . Ultracold plasmas not 
only allow for experimental tests of present theories but also for studies of the 
ionic correlation dynamics in expanding systems with a steadily changing local 
equilibrium state, as described by the hydrodynamical evolution equations. 
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The modified time evolution of the ionic temperature leads to a modified ex- 
pansion dynamics of the plasma through the change in the ionic pressure, 
which contains a correlation part in addition to the thermal part. The devel- 
opment of ionic correlations also reduces of the total electrostatic interaction 
energy compared to a fully uncorrelated plasma. Together, these two effects 
result in the term U\\/?> on the rhs of f!39bj) . which constitutes an effective 
negative acceleration slowing down the plasma expansion. This term can be 
regarded as a time dependent "external" force, giving rise to an additional ef- 
fective potential in which the ions move. This potential, however, is not static, 
but changes during the expansion of the plasma, and leads (due to conserva- 
tion of total energy) to an additional heating of the ions (second term on the 
rhs of PS) ). 



Eqs. f l39al) - fl39cp . f!31dl) and (13 lei) form the final set of equations for the macro- 
scopic parameters of the plasma state. They are quickly and easily solved nu- 
merically. Since the numerical effort is independent of the number of particles 
involved, the hydrodynamic approach provides a convenient method for simu- 
lations of large plasma clouds over long, i.e., millisecond, timescales. It can be 
used to quickly gain insight into the plasma dynamics by efficiently scanning a 
broad range of initial-state parameters, and it is able to simulate plasmas that 
are too large to be treated effectively with molecular dynamics methods. More- 
over, and maybe even more importantly, the kinetic model provides physical 
insight by reducing the plasma dynamics to a few macroscopic parameters. 
However, the accuracy of these approximate macroscopic treatments, partic- 
ularly when they describe dynamics influenced by ion correlation, must be 
assessed carefully. This can be done to some extent by comparison with ex- 
periment, but more rigorously by developing numerically accurate microscopic 
approaches relaxing most of the approximations invoked in the macroscopic 
approach. 



3.2 Microscopic Approaches 



Microscopic approaches are by definition more accurate. They are of great 
value for justifying macroscopic descriptions, and of course are the only al- 
ternatives for plasma conditions that preclude macroscopic descriptions. This 
refers to detailed information about the very early system evolution which is 
not accessible to hydrodynamics, and the effects of strong particle correlations, 
which are of particular interest in the context of ultracold plasmas. 
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3.2.1 Molecular Dynamics Simulations 

Molecular dynamics simulations of Coulomb-type interacting particles have 
been employed for investigating a variety of different problems in plasma 
physics and are a corner stone of astrophysical studies of galaxy and glob- 
ular clusters |105j . Conceptually, the problem seems rather simple as it only 
involves the solution of a set of coupled differential equations describing the 
classical motion of charged particles. However, there are several technical chal- 
lenges, which in the case of ultracold neutral plasmas arise from the open 
boundary of the expanding system, the attractive nature of the electron-ion 
interaction and the small mass ratio of the electrons and the ions. 

A general difficulty for simulations of neutral plasmas arises from the bound 
state dynamics of atomic electrons, whose orbital timescale is many orders of 
magnitude smaller than the plasma evolution time. In particular for highly 
eccentric orbits conventional integrators, such as higher order Runge-Kutta 
methods, do not provide a stable solution, as they lead to a significant drift 
of the total energy. More appropriate propagation schemes such as symplectic 
integrators without energy drift and originally designed for the integration 
of gravitational systems, have clear advantages, but also several draw-backs 
and the method of choice depends on the particular problem to be studied. 
For simulating ultracold plasmas higher order symplectic integrators [106j . 
predictor-corrector schemes [JT] arid timescaling methods combined with reg- 
ularization procedures |lU7j can yield stable and reliable results. 

Although a single electron and ion pair can be propagated very stably, many- 
body simulations have to fight additional complications. Electrons are bound 
to different ions with different binding energies and hence their orbital times 
can be very different. Moreover, these orbital times are generally much smaller 
than the timescale on which the heavier ions move. With a typical ion timescale 
on the order of microseconds, one has to cover timescales from femtoseconds to 
microseconds in order to describe the initial relaxation of the ionic component. 
Mazevet et al. [38] and Kuzmin and O'Neil [JT] have used a small ion mass and 
scaled the obtained timescales to realistic experimental situations to obtain 
a quantitative estimate of the experimental results. Such a scaling, however, 
does not apply to all aspects of the plasma dynamics, since e.g., the timescales 
for electron-ion equilibration and ion-ion equilibration scale with a different 
power of the ion mass. 

Several methods to numerically tackle very different timescales have been dis- 
cussed in the literature (see, e.g., [1U8] and |1U5|I109| for an overview and 
further references). The common idea of such approaches is to introduce in- 
dividual timesteps for the different particles. The criteria determining the 
timestep and the procedures for the synchronization of the different particles 
can however be quite different and have to be adapted to the specific prob- 
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lem under consideration. Kuzmin and O'Neil [^T|37] used a modified version 
of the method originally developed by Aarseth [1 lOj to perform a full scale 
simulation of ultracold neutral plasmas with a realistic xenon ion mass and a 
total number of 8192 particles. 

The number of particles N is generally limited by the number of force evalu- 
ations that have to be carried out for every timestep. Since for every particle 
one has to sum up N — 1 forces of the surrounding charges the numerical effort 
quickly increases as iV • (N — 1) ~ iV 2 with increasing number of particles. 
Simulations of large, extended plasmas typically consider a smaller part of the 
system and employ periodic boundary conditions to mimic the properties of 
the larger system |109j . For ultracold plasmas such an approach is not very 
appropriate as it neglects an important part of the dynamics - the expansion 
of the plasma. The A r2 -scaling of the numerical effort can, however, be re- 
duced by so-called hierarchical tree-methods |lllj or fast multipole methods 
[112] . originally developed to simulate globular clusters. Such approaches still 
account for the exact forces of nearby charges but approximate the force of 
distant particles by the average force up to a certain multipole order of the 
corresponding charge distribution. This reduces the scaling of the numerical 
effort with the number of particles to ~ iV In iV or even ~ N. 



3.2.2 Particle- In- Cell Approaches 

An alternative and conceptually different method to reduce the numerical ef- 
fort connected with the summation of the inter-particle forces is to neglect cor- 
relations between the charges but retain their meanfield interaction due to the 
long-range part of the Coulomb potential. The so-called particle-in-cell simu- 
lations (PIC) |113j have become a standard method to treat various plasma 
physics problems, such as the collective dynamics of fusion devices jl 14(1115] 
or of cold plasmas in gas discharges |116|I117|I118] . i.e., weakly coupled plas- 
mas. Clearly, for ultracold neutral plasmas this method seems questionable, 
as we are interested in the regime of strong coupling. However, the approach 
has proven to yield a reliable and efficient description of certain aspects of 
the plasma dynamics as long as the plasma expansion is mainly driven by the 
thermal pressure of the electrons [39J. 

The striking advantage of the PIC method is the possibility to treat very large 
systems, by replacing a collection of real particles by a single super-particle, 
with an increased charge and mass, retaining the charge-to-mass ratio of the 
physical particle. A single propagation timestep of such an ensemble of super- 
particles consists of the following sub-steps: 

(1) At the beginning of each propagation timestep the positions of these 
super-particles are used to determine the charge density of the system, 
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represented on an appropriately chosen grid. To reduce numerical noise 
due the finite particle number the charges are represented by a shape 
function. Most common shape functions are rectangles (particle in cell), 
tent functions (cloud in cell) or Gaussians centered around the actual 
particle position. Further problems arise for particle weighting in curved 
coordinates, as discussed in |119j . 

(2) Having obtained the charge density one evaluates the corresponding static 
electric field from the Poisson equation (Eq. [T9lh If the system exhibits 
additional spatial symmetries this is easily done by employing Gauss' 
theorem. In the general case the electric mean-field is most efficiently 
evaluated within an FFT procedure. 

(3) The calculated electric field is subsequently transformed to the actual 
particle positions taking into account the chosen shape function, which 
finally allows advancement of all particle positions and velocities accord- 
ing to the classical equations of motion. 

The positions and velocities of the test-particles have no physical meaning. The 
relevant quantities are the phase space distributions of the charges determined 
by the ensemble of test-particles. It is easy to see that the resulting dynamics 
provides a numerical solution of the Vlasov equation (fT8"j) . 

Hence, the PIC treatment as described so far does not account for any colli- 
sional processes, such as formation of Rydberg atoms and subsequent electron- 
atom scattering. Such collisions can however be incorporated into the PIC- 
treatment on the basis of the collision integrals introduced in Eqs. ( !2T|) - 
(1291) . Within such so-called particle- in-cell-Monte-Carlo-collision simulations 
(PICMCC) [86] the collision terms in Eq. ([TBI are evaluated by random real- 
izations of collision events, with a probability distribution determined by the 
respective collision integral. 

For the collision integrals given in Eqs. f[2"8l and f[2U|) describing inelastic col- 
lisions and Rydberg atom formation, this stochastic evaluation procedure is 
particularly simple. First, one introduces an additional species by providing 
each ion with an additional label n, which characterizes its internal state. 
If the particle is an ion, then n = 0, and n > for atoms, where n de- 
notes the principle quantum number of the atom. During each timestep of 
length At, an ion at position r is transformed into an atom with probabil- 
ity P t br = At J2n Ktbr(n, p e (i"i), T e ) . In this case, the corresponding princi- 
pal quantum number is sampled from the probability distribution P t b r (n) = 
AtK thr (n, p e (r{) ,T e ) / P thr . At the same time the electron number is decreased 
by unity. Analogously, for bound states, an electron-atom collision occurs 
with the probability P ea = At(K ion (n, p e (r a ), T e ) + T, PJ t n K hh (n,p, p e (r a ), T e )). 
With a probability Pi on = AtK ion (n, p c (r a ), T e )/P ea , this leads to ionization 
of the atom. Otherwise, the collision results in a bound-bound transition 
to a state with principal quantum number p with probability Pbb(^,p) = 
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At K hh (n, p, p e ( r a ) , T e ) /P ea . 

We finally note, that the tractable collision integrals are not limited to such 
a simple case, but can also be nonlinear in the distribution functions. In par- 
ticular, elastic electron-ion and ion-ion collisions can be described in a similar 
way p3M] . 



3.2.3 Hybrid Methods 

So far we have discussed two complementary approaches which either provide 
the full microscopic information about the early plasma evolution or allow 
to follow the long-time dynamics on the cost of neglecting possible strong 
correlations between the plasma particles. 

While the ionic component may indeed form a strongly correlated state, it 
turns out that the electrons remain weakly coupled during the entire plasma 
evolution (see Section [4.51 and 14.61) . In [34"|39] this different behavior of the 
two particle species has been used to develop a hybrid-molecular-dynamics 
approach. It treats the electrons and ions on different levels of approximations, 
combining the advantages of a mean field treatment (electron dynamics) and 
molecular dynamics simulations (ion motion). Collisions between electrons 
and the heavy particles (atoms and ions) are accounted for within the Monte- 
Carlo procedure described in the previous subsection. By eliminating the short 
dynamical timescale of the bound electron motion this method significantly 
reduces the numerical workload, while simultaneously enabling a virtually 
exact numerical simulation of the ion dynamics fully incorporating ion-ion 
correlations. 

The numerical effort may be further reduced within an adiabatic approxima- 
tion for the electrons. As discussed in Section 13.1.11 the small electron-to-ion 
mass ratio permits an adiabatic approximation, where instantaneous (on the 
timescale of the ion dynamics) equilibration of the electrons is assumed. Keep- 
ing this approximation in a (hybrid) MD approach means that much larger 
timesteps may be chosen, since the electronic dynamics does not need to be 
resolved. In this way, an MD simulation of the ion motion over the whole 
timescale of the experiments becomes feasible. According to Eq. ([20]) . the elec- 
tronic distribution function is then given by a quasistatic distribution at 
all times. 

The exact form of f^ qs \ however, merits some discussion. The assumption of 
an ideal equilibrium, i.e. a Maxwell-Boltzmann distribution 




(47) 
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leads to fundamental problems. Since the mean-field potential ip has a finite 
depth, the corresponding Maxwell-Boltzmann distribution approaches a finite 
non-zero value even for r — > oo and hence is not normalizable. In |33j, this 
problem was circumvented by introducing a spatially dependent electron tem- 
perature which is constant within a finite region around the plasma center, 
but decreases quadratically for large distances. An analogous problem was dis- 
cussed long ago in an astrophysical context, in connection with the formation 
of star clusters, where the assumption of a Maxwellian velocity distribution 
implies an infinite mass of the star cluster in complete analogy to the present 
case. As discussed by Chandrasekhar [120], such a velocity distribution con- 
tains particles (stars) with arbitrarily high velocities which are not bound to 
the cluster potential. The rate with which these stars evaporate from the clus- 
ter has also been considered in [120j . and later in [121] . If the evaporation is 
much slower than the relaxation of the stars remaining in the cluster, a qua- 
sistationary distribution forms. Comparison with numerical simulations shows 
that this distribution is well described by a truncated Maxwell-Boltzmann dis- 
tribution 
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in the center of the cluster. The escape velocity w e = 
by the depth of the potential well, S(p = <p(r — > oo 
this may be generalized to arbitrary positions, 
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If the potential is a monotonic function of the distance from the cluster center, 
as in the case of gravitational forces, the space-dependent escape velocity is 
given by w e - 



2e[(f(oo) — (p(r c )]/m e for a spherical symmetric system, and 
one obtains the so-called Michie-King distribution [123|IT2 2j. In the present 
case of a neutral plasma, the potential may be non-monotonic due to the 
presence of charges with different sign. Hence, the escape velocity has to be 
defined more generally as 
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With W e = m e w\j (2k B T e ), the electronic density is then given by 
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which approaches a Maxwell-Boltzmann distribution in the center of the plasma 
for sufficiently low temperatures (i.e. W e 3> 1) but decays significantly faster 
towards the edge of the plasma. 
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The Poisson equation, particle number conservation, energy conservation, and 
distribution Eq. (|5T|) constitute a closed set of equations 



A<p = Ane (p e — Pi) (52a) 

Peocexpf-^) [ Wc e~ x x 3/2 dx (52b) 
\k B l e J Jo 

N e = J p e dr e (52c) 

E ges = -7r vlf e dv e d\ e - \ p c ipdr e + E { + E a = const. , (52d) 



that uniquely determines the electronic density, temperature, and mean-field 
potential <p at all times. (In Eq. (152dl) . is the sum of potential and kinetic 
energies of all ions, and E a is the sum of the electronic binding energies and 
kinetic energies of all atoms.) Having calculated the electron density profile at 
time t, the knowledge of the corresponding electron mean field force is used to 
advance the ion positions and velocities under the influence of the electronic 
mean field and the full ion-ion interaction according to 

m{Y- hj = F ei + F ii ( r i,j) r i,k) , (53) 

where ry is the position of the jth ion, F e i is the force due to the electronic 
mean-field, 

F ei = -e^ ; Av? c = 47rep e , (54) 

and 

F ii (r i „r i , fc )= ri ' J ~ ri 'V (55) 

In this propagation scheme it is still numerically quite costly to calculate the 
inter-ionic forces. However, as demonstrated in [30] , with usage of hierarchical 
tree-methods for the direct propagation of the ions, this hybrid-MD approach 
allows the description of large systems, with ~ 10 6 ions over simulation times 
of several tens of microseconds. 



4 Physical Processes in Ultracold Neutral Plasmas 

The previous section presented the experimental and theoretical tools to gen- 
erate, detect, and understand ultracold neutral plasmas. In this section we will 
use these tools to discuss physical processes. We proceed along the timeline 
of the dynamics (Fig. H]), i.e., we start in Section 14.11 with electron equili- 
bration after photoionization. Then we will cover the establishment of local 
thermal equilibrium for the ions in Section 14.21 Before we come to the plasma 
expansion (Section 14.41) into the surrounding vacuum, we introduce collective 



33 



&0 



a 



a 
o 



















\ f 






\ t 1 






V_i_y 







Position 






Fig. 9. Schematic of the potential energy seen by a test electron when enough atoms 
are photoionized to result in trapping of electrons. Photoionization occurs at to = 
and the sample is neutral everywhere. Because of the kinetic energy imparted by 
the laser, some electrons leave and a charge imbalance develops. At t\ ~ 10 ns the 
resulting potential well equals the initial kinetic energy, trapping the remaining 
electrons. As electrons in the well thermalize, evaporation occurs. The well depth 
increases and the electrons cool slightly. By t<i « 1 //s evaporation essentially stops. 
The dashed line indicates the average kinetic energy of the electrons. Reused with 
permission from [2]. Copyright 1999, American Physical Society. 

excitations of the plasma electrons, which may be used to probe the plasma 
expansion dynamics (Section |4.3p . 



Two more subsections are devoted to important topics. First, a detailed analy- 
sis of the initial heating of the electrons, which affects the expansion dynamics 
is given in Section 14.51 Secondly, Section 14.61 will discuss in some detail the 
evolution of the electronic and ionic Coulomb coupling parameters, which are 
important since they indicate the degree of correlations that may be realized 
in a plasma. 



4-1 Initial Electron Equilibration 



4-1.1 Trapping of Electrons 

The first experimental study of ultracold neutral plasmas [2] explored the cre- 
ation of the plasma and trapping of the electrons by the ionic background 
charge. (See Fig. El) Immediately after photoionization, the charge distribu- 
tion is neutral everywhere. Due to the kinetic energy of the electrons, the 
electron cloud expands on the timescale of the inverse electron plasma fre- 
quency r e = u~l = ^m e e / p e e 2 ~ 1 ns, where p e is the electron density, and 
cjp e is the electron plasma oscillation frequency. On this timescale the ions are 
essentially immobile. The resulting local charge imbalance creates a Coulomb 
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Fig. 10. Electron signals recorded for four different pulse energies of the photoion- 
izing laser, i.e., different densities of charged particles (10 5 — 10 7 cm~ 3 ). The upper- 
most curve corresponds to the lowest energy. The photoionization occurs at t = 0. 
The initial kinetic energy of the electrons is AE/ks = 0.6 K. The data shown is 
an average over 20 cycles of the experiment. Also shown is the magnitude of the 
applied electric field. Reused with permission from [2]. Copyright 1999, American 
Physical Society. 



potential energy well that traps all but a small fraction (< 5%) of the electrons. 
Smaller clouds, more ions and electrons, and lower initial electron kinetic en- 
ergy all lead to a larger fraction of electrons trapped by the ions. Simple PICC 
simulations [2] show that electrons escape mostly from the edges of the spatial 
distribution. 

This description implies that there is a threshold for trapping electrons. If 
the well never becomes deeper than the initial kinetic energy, all the electrons 
escape. This is demonstrated in Fig. []j]J and the diagnostic used to obtain 
this data was charged particle detection of electrons after they have left the 
plasma, as shown in Fig. [5j In each cycle of the experiment, the atoms are 
first laser cooled and an electric field of approximately 5mV/cm is applied. 
Subsequently, the atoms are photoionized, and after about 500 ns of time of 
flight a pulse of electrons arrives at the detector (see the first peak in the curves 
of Fig. [T0l) . A few microseconds later, the electric field is linearly increased. If 
electrons have been trapped by the ionic background potential, they will be 
liberated now to produce the second peak in the curves of Fig. [TU1 If the laser 
intensity is so low that the well never becomes deeper than the initial kinetic 
energy, no electrons are trapped and the second peak is missing (top trace). In 
addition, when electrons are trapped, the first peak develops a tail. Trapped 
electrons thermalize within 10 — 100 ns [60], and as charges are promoted to 
energies above the trap depth, they leave the well. On the time scale of this 
experiment, the ions are essentially stationary. 
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Fig. 11. N* , threshold number of ions required for trapping electrons for a Gaussian 
plasma with a = 1 mm. 

One can quantify the threshold for the trapping effect. The last electron to 
leave the plasma must climb out of the potential formed by the ions, with 
no neutralizing electrons. At threshold, the depth of this well equals E e . The 
density profile of the ions follows that of the neutral atoms Eq. (j2J), where 
the maximum density in the center is given by p\ = A^/(27r(T 2 ) 3//2 with the 
total number of ions, iVj. The electrostatic potential produced by this charge 
distribution is 



U(r,N$ = - 



47r£n|r — r 



T pi(r'W = —TV; erf 



(56) 



where erf (re) is the error function. The depth of the potential well is given by 
its value at the center r = 0, 
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When — £7(0, N{) equals the electron kinetic energy, E e , then the last electron 
to leave the well is trapped. In other words, the minimum number of ions 
required to trap any electrons, defined as N*, is found from E e = —U(0, N*). 
Figure [TT] shows typical values of iV* . 

The threshold condition iV; = iV* is mathematically equivalent to Ad = a, 
where A^> is the Debye screening length as defined in Eq. f|45|) with T e = 
|£' e //cB- Note, that the relation a > X D is the standard condition for a sys- 
tem to be considered a plasma [124] . namely that the system size (a) has to 
be larger than the characteristic length scale (Ad) on which local charge im- 
balances are screened. In the present context, we can also interpret A^> as the 
displacement of electrons from their equilibrium positions when their energy in 
the local internal electric field in the plasma equals their kinetic energy [125j . 
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Fig. 12. (left) The fraction of electrons trapped is plotted versus the number of 
photoions created. Each curve corresponds to a different (green) laser frequency. The 
corresponding initial energies of the electrons are displayed in the legend, (right) 
Same as (left) but the number of ions is scaled by N*, see text. There is a scale 
uncertainty of about 10% in determining the fraction of electrons trapped. Reused 
with permission from [2j. Copyright 1999, American Physical Society. 

If Xd > cr, the electrons are free to escape to infinity. If A^> < a, electrons are 
trapped by the ion cloud. In the center of a typical ultracold neutral plasma, 
X D = 1 — 10 yum. Typical plasma sizes are on the order of some hundred jttm, 
such that the requirement for the existence of a plasma state is well fulfilled 
under experimental conditions. 

Fig. [12] [2 J illustrates the trapping effect and the role of N*. We see the number 
of ions required to begin trapping electrons grows with increasing E e (left 
panel), i.e., the threshold N* increases. When the number of ions is scaled 
by N*, the fraction of electrons trapped follows a universal curve, clearly 
exhibiting the threshold N±/N* = 1 (right panel). For iVj < N* no electrons 
are trapped. This figure also shows that for N/N* ^> 1, the non- neutrality is 
small with the fraction of trapped electrons almost reaching unity. 

The fraction of trapped electrons is surprisingly well described by the simple 
formula [3711] 

N e _ y/NjN* - 1 



An analytic derivation of this dependence has not yet been found. 



(5* 



A theoretical discussion assuming a quasi-neutral plasma gives a better un- 
derstanding of the degree and nature of non-neutrality in the plasma. First 
note that the quasi-neutrality condition does not imply an exact equality of 
the electron and ion densities and hence a vanishing space charge potential ac- 
cording to the Poisson equation (fl9l) . Instead, the Poisson equation is replaced 
by the electron Vlasov equation which in adiabatic approximation yields for 



37 




Fig. 13. (a) Electron velocity distribution at t = 40ns (dash-dotted line), t = 80ns 
(dotted line) and t = 160 ns (dashed line) compared to a Maxwell distribution (solid 
line) . The distribution has been sampled from a shell 1 /4 out from the center of the 
plasma, (b) Time evolution of the electron temperature at the plasma center (dotted 
line), 1/4 out from the center (dashed line), 2/4 out from the center (dash-dotted 
line) and 3/4 out from the center (dash-dot-dot-dot line). The solid line shows 
the average temperature of the innermost 80% of the electrons. The initial state 
parameters are p e = 10 9 cm -3 , T e = 100 K and a = 196 pm. Reused with permission 
from [33]. Copyright 2003, American Institute of Physics. 



the space charge potential (Eq. fl2T1) ) 

e<p(r) = -k B T c In Pi(r) = h B T e ^ . (59) 
Pi{r = 0) 2 <y l 

Now we can use this expression for a backward calculation of the charge sep- 
aration from the Possion equation, which gives a spatially constant excess 
charge of 

k\i 

pi - p c = 3 (60) 
e z a z 

Introducing the local Debye length, which is a function of density and thus 
radius, we find from Eq. (I60l . that the relative charge difference, i.e., the local 
degree of non-neutrality 

Pi ~ Pc ^lij) t f\1 \ 

= — — , 61 

Pe cr 2 

is solely determined by the ratio of the Debye length and the system size — just 
like the global non-neutrality (Eq. (158!) ). The fractional non-neutrality is small 
as long as \d < a, typically less than 1% in the center of the plasma, where 
the density is high. It becomes large near the edge, however, violating the 
quasi-neutrality assumption, which may lead to deviations from the Gaussian 
density profile (see Section l4.4.4l for a more detailed discussion). For the major 
parts of the plasma volume A^(r) < a 2 , justifying the statement that the 
plasma is essentially in a quasineutral state. 

Non-neutrality is also related to a spatially dependent electron temperature. 
When using the definition of the Debye length Ad in the above argumentat 
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we have implicitly assumed the existence of an electronic temperature. While 
this is certainly a good approximation, a global thermal electron velocity dis- 
tribution is just developing during the initial phase of the plasma relaxation. 
PICMCC simulations (Section 13.2.21) . which yield very good agreement with 
the measured evaporating electron flux [33], show that the velocity distribu- 
tion in the central plasma region quickly reaches its Maxwellian form within 
a few hundred 100 ns (see Fig. [T3a). As shown in Fig. fl~3b it takes consider- 
ably longer for a global equilibrium to develop. The center of the plasma is 
seen to heat up, while the outer region cools down. The total electron kinetic 
energy is, however, almost unaffected by these local temperature imbalances. 
In reality, the electron gas heats up due to the formation of Rydberg atoms 
[66,42,33,39] and an initial build-up of spatial electron correlations. The latter 
process takes place simultaneously with the initial electron evaporation, but 
is not observable in Fig. [L_31 due to the inability of the PIC treatment to de- 
scribe the correlated dynamics of the electrons. A more detailed discussion of 
the various electron heating mechanisms will be given in Section 14.51 

4-1-2 Spontaneous Ionization of a Dense Rydberg Gas 

Compared to photoionization of a trapped cloud of ground state atoms, the 
dynamics of plasma formation is more complicated when it is created by spon- 
taneous ionization in a dense cloud of ultracold Rydberg atoms. Reference [56] 
showed that some Rydberg atoms are initially ionized by black-body radia- 
tion, collisions with hot background atoms, or Penning collisions with other 
Rydberg atoms. Attractive dipole-dipole and van der Waals interactions that 
pull Rydberg atoms together and enhance the Penning ionization rate are 
particularly important for high principal quantum number (n > 50) [126J. 
Resulting electrons escape the cloud in the first few microseconds until a few 
thousand ions build up to form a Coulomb well that traps subsequently pro- 
duced electrons. This is reminiscent of the model presented in Fig. At this 
point, the electrons seed an ionization avalanche. As much as 2/3 of the ini- 
tial Rydberg population is ionized, and the remaining fraction is collisionally 
de-excited to lower Rydberg levels, satisfying energy conservation [55]. Dur- 
ing the avalanche, /-changing electron- Rydberg collisions, n-changing electron- 
Rydberg collisions, electron-Rydberg ionizing collisions, and Penning ionizing 
Rydberg- Rydberg collisions all play a role [55|127j . By exciting Rydberg atoms 
with a relatively low density in a pre-existing ultracold plasma, it was shown 
in [68] that very tightly bound Rydberg atoms (n ~ 20) could be ionized by 
electron-Rydberg collisions in these systems. The energetics of electrons in 
plasmas created in this way, such as the equilibrium temperature after the 
first few tens of nanoseconds, is a question that has not been well-explored. 
Theory [33] predicts no drastic change as the excitation laser crosses the ion- 
ization threshold and the experiment changes from direct photoionization to 
initial excitation of Rydberg atoms. Recent experimental work supports this 
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Fig. 14. Time evolution of the effective ion temperature for a Sr plasma with initial 
peak density po(0) = 2 ■ 10 9 cm -3 and electron temperature T e (0) = 38 K. The 
results of a hybrid-MD calculation (solid line) are compared to experimental data 
extracted from an absorption imaging measurement (dots). Reused with permission 
from [30J. Copyright 2005, American Physical Society. 

prediction [70] . 

4-2 Initial Ion Equilibration 

The establishment of LTE for the ions occurs on a longer timescale than it 
does for the electrons, and it is conveniently studied experimentally. In addi- 
tion, several inelastic processes (most notably three-body recombination) that 
heat the electrons do not affect the ions. In this unusual regime, interesting 
phenomena are expected to be observable. 

Theoretically, the ion dynamics can be studied within the Hybrid Molecular 
Dynamics approach outlined in Section 13.2.31 In this approach, the ions are 
propagated individually, so that full information about their positions, veloc- 
ities, as well as spatial correlations is available. Experimentally, on the other 
hand, the optical techniques described in Section 12.2.21 provide an excellent 
tool for following the ion dynamics. 

4-2.1 Disorder-induced heating 

Figure [TH shows that the experimental evolution of the ion temperature of an 
ultracold plasma and the theoretical curve obtained from a hybrid-MD simula- 
tion are in good quantitative agreement over a time relevant for ion thermaliza- 
tion |30j. The most striking observation from Fig. [13] is the rapid heating of the 
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plasma ions by several orders of magnitude, which is the so-called "disorder- 
induced heating" or "correlation heating." It originates from the fact that the 
system is created in an undercorrelated state far from thermodynamic equi- 
librium. Immediately after photoionization, the ions in an ultracold neutral 
plasma have very little kinetic energy. But they are spatially uncorrelated, 
and there is significant excess potential energy in the system compared to the 
equilibrium state in which spatial correlations keep ions from being close to 
each other |128lfTT] . During early evolution of the plasma, correlations develop 
and potential energy is converted into (thermal) kinetic energy [4"0"|103|I104] 
leading to equilibration. However, during the disorder-induced heating phase, 
the plasma is not in LTE, and the notion of an ion "temperature" should be 
interpreted as a measure of the average ion kinetic energy. 

Within the framework of the hydrodynamical model of Section 13.1.31 the 
above interpretation of disorder-induced heating is immediately apparent in 
Eq. fl39cl) : As long as expansion of the plasma is still negligible (7 = 0), the 
change in ion temperature is directly determined by the build-up of ion-ion 
correlations (last term on the right-hand side of f!39cj) ) . This process has been 
discussed in theoretical papers in many contexts. Early interest was generated 
by non-equilibrium plasmas created by fast-pulse laser irradiation of solid 
targets [129,130 : 131,80. 8"1"|132| . Experimental results were lacking, however, 
because of the fast timescales involved and limited possibilities for diagnos- 
tics. Ultracold neutral plasmas - much better suited for a detailed study of the 
process - reinvigorated interest in the mechanism of disorder-induced heating 
[3TllM0lfT03llU4lll33j . 

Qualitatively, one expects the timescale for the ion heating to be given by 
the inverse of the ionic plasma frequency, which is the typical timescale 
on which spatial correlations develop in a plasma. Physically, co~l is the time 
for an ion to move one interparticle spacing when accelerated by a typical 
Coulomb force of e 2 /47re a 2 , where a = (47rpi/3) -1 / 3 is the Wigner-Seitz ra- 
dius. For the typical plasma of Fig. [TH a;" I « 270 ns in good agreement 
with the rise-time of the temperature observed in the figure. This timescale 
of ~ 10 2 ns is several orders of magnitude larger than in laser-generated high- 
density plasmas as mentioned above, and it is easily resolved experimentally. 
One would guess the ion temperature after equilibration would be on the order 
of the Coulomb interaction energy between neighboring ions, 



e 2 1 



47T£o a 



U a , (62) 



which is a few kelvin. A quantitative analysis |40j, assuming complete initial 
disorder and incorporating the screening effects of the electrons, predicts from 
a simple energy conservation argument an equilibrium ion temperature of 
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Fig. 15. The effective ion temperature, T i e fj (Eq. (fl2j) ). versus time after photoion- 
ization for initial electron temperature of T e = 2E e /3kB = 38 ± 6K and various 
plasma densities, (a) The data is plotted on absolute temperature and time scales, 
(b) The time is scaled by the inverse of the average plasma period, and lj )e g is 
scaled by Tq (Eq. [62]). Reused with permission from [28J. Copyright 2004, Ameri- 
can Physical Society. 
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Here, k = a/Xo is the Debye screening constant, and U is the excess potential 
energy per ion. U has been studied with molecular dynamics simulations |134] 
for a homogeneous system of particles interacting through a Yukawa poten- 
tial, <f)(r) = U a (a/r) exp(— r/X D ), which describes ions in the background of 
weakly coupled electrons [135j . To obtain 7] for given p x and T e , one needs to 
selfconsistently solve Eq. (1631) with an analytic expression for U [102J since U 
depends itself on 2]. 



The arguments given above imply that both the timescale and the energy scale 
of the disorder-induced heating depend on the plasma density only (apart from 
a small dependence of the final temperature on the screening k, which will be 
discussed below). Fig. [T5l confirms this expectation. If the measured effective 
ion temperature T i eff is scaled with T c (Eq. (1621 ) and time with the inverse 
plasma frequency u p ^ the measured curves for different densities coincide quite 
well (Fig. [ToT b)). Small differences in the curves are due to electron screening 
of the ion-ion interaction. The effect of electron is demonstrated in Fig. [T6| 
where Ti jC fj obtained from absorption spectra (Eq. (|12p ) is shown for different 
initial electron energies but the same ion density distribution: Smaller initial 
electron energy (temperature) implies more screening and a lower equilibrium 
ion temperature. The latter is directly related to the initial excess potential 
energy (Eq. (|63|) ). which is lowered if screening through electrons is increased. 
It was found that the measured ion temperature just after the plasma reaches 
LTE agrees well with the theoretical value T iiave according to Eq. (|63l) and 
averaged over the ion distribution 
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Fig. 16. Effective ion temperature versus time after photoionization for var- 
ious initial electron kinetic energies. The peak ion density for all data is 
poi = (1.4 ± 0.5) x 10 10 cm~ 3 . Reused with permission from [28J. Copyright 2004, 
American Physical Society. 

As emphasized in |133j . it is not surprising that that Eq. fl63|) accurately pre- 
dicts the ion temperature for equilibrating ultracold neutral plasmas, because 
it basically expresses energy conservation. However, it is worth pointing out 
that for data in Ref. [28] with the lowest fit T e and highest poi, the peak value 
of k in the plasma is 0.7. This corresponds to three electrons per Debye sphere 
(k~ 3 = p c 47rA 3 :) /3). One might not necessarily expect Eq. (|63|) to be accurate 
in this regime because it assumes a Yukawa potential for ion-ion interactions, 
which is based on Debye screening and is normally derived for k -3 ^> 1. 



4-2.2 Kinetic energy oscillations 

A close inspection of Figs. [15] and [16] reveals that at the end of the disorder- 
induced heating phase, the ion temperature overshoots its equilibrium value 
before settling to it. This phenomenon is even more evident in Fig. [T7|, where 
Ti >e ff is shown for an inner and outer region of the plasma image (p = \/x 2 + y' 2 < 
0.9 cr and p > 1.48 a, respectively). Each selected annular region contains 1/3 
of the ions and probes a region with significantly less variation in density than 
in the entire plasma. The region with lower density has lower ion temperature, 
as expected from Eq. fl6"3"|) . but the oscillation is the most striking observation. 
The oscillation was also observed in a calcium plasma with fluorescence mea- 
surement [13] (Fig. [T8l) . The contrast of the oscillations is strong in this case 
because the small volume of the plasma excited by the fluorescence excitation 
laser significantly reduced the variation of the sampled density. Numerical 
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Fig. 17. Effective ion temperature obtained from different selected regions of a 
plasma cloud with p oi = (4 ± 2) x 10 9 cm" 3 and initial T e = 2E c /3k B = 38 ± 6K. 
Oscillations are clearly visible at the end of the disorder-induced heating process. 
Reused with permission from [28J. Copyright 2004, American Physical Society. 




°0 0.25 0.5 0.75 1 

time (jis) 

Fig. 18. The mean z-component of the ion velocity in a calcium plasma as a function 
of time after photoionization in a region defined by the excitation laser beam. The 
velocity was extracted using Eq. (|17p . and the absolute calibration was obtained by 
assuming the ions had an initial velocity distribution equal to that of the laser-cooled 
atoms. The rapid increase of the velocity at early times is due to disorder-induced 
heating. The subsequent slower increase reflects the expansion of the plasma, which 
will be discussed in Section 14.41 The solid line is a fit of the velocity to a model for 
the plasma expansion plus an offset representing the random thermal ion velocity. 
The inset shows the ion velocity divided by the fit, which highlights the oscillation 
of the kinetic energy. Reused with permission from [13]. Copyright 2005, American 
Physical Society. 

simulations clearly show the oscillations and also provide their temperature 
and spatial dependence (Fig. [191) . 

These oscillations display universal relaxation dynamics of a strongly coupled 
Coulomb system. More precisely, they reflect the fact that the two-particle 
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Fig. 19. Spatio-temporal evolution of the ion temperature, obtained from HMD-sim- 
ulations for 40000 ions and electrons with an initial Coulomb coupling parameter 
of r e (0) = 0.07. Reused with permission from [29]. Copyright 2005, Institute of 
Physics. 

distribution function (i.e. spatial correlations) and the one-particle distribu- 
tion function (i.e. temperature) relax on the same timescale (l/u Pi i) due to the 
strong long-range interactions. Normally, correlations relax much more quickly, 
and the separation of these time scales is known as the Bogoliubov assumption 
[136J. Studies with ultracold neutral plasmas represent the first experimental 
observation of kinetic energy oscillations, but it has been the subject of in- 
tense theoretical study through analytic calculations [130j and simulations 
[129j[ 131 80 : 81~yi~32l)133f30j of equilibrating strongly coupled plasmas. Calcula- 
tions |129|137j also show oscillations in the velocity autocorrelation function 
in equilibrium systems. 

Intuitively, one may relate this phenomenon to the motion of each ion in its 
local potential energy well. The discussion of disorder-induced heating in the 
previous section qualitatively ascribed the increase in ion kinetic energy to 
each ion moving towards the minimum of its well. The time scale for this 
motion is uj~1. However, the ion will not suddenly stop at the bottom of its 
potential; it will overshoot and climb up the hill again. It is reasonable to 
expect harmonic motion to persist for some time, as kinetic and potential 
energy are interchanged. This suggests a kinetic energy oscillation at 2u; Pj i, 
which is in line with the observations. Also consistent with this description, 
the oscillation period observed experimentally in Fig. [T7] is longer in the outer 
region where the average density is lower and hence the plasma frequency is 
smaller. This also explains why averaging over the entire cloud such as in Fig. 
[141 obscures the oscillation: the motion dephases because of the variation in 
o; p> i. In a disordered system corresponding to the ultracold plasma scenario, 
there is probably no collective or long-range coherence to the motion, so in 
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spite of the fact that u Pj i sets the timescale of the oscillation, one must be 
cautious about describing the motion as an ion plasma oscillation. 

The numerical hybrid-MD approach permits the definition of a local temper- 
ature. Hence, it provides better temporal and spatial resolution of T\ than the 
experimental optical probes. As was shown in [73 29], the ion temperature 
shows clear oscillations in space (i.e. as a function of the distance from the 
center of the plasma) and, in fact, a wavelike oscillation pattern emanating 
from the plasma center (Fig. [T91) . 

From the experimental data, it is not straightforward to comment on the 
damping of the observed kinetic energy oscillations because even in Fig. [171 
which provides an annularly-resolved, more-local probe of the system, the 
analysis still averages over the density variation along the z-axis of the plasma. 
This introduces dephasing because the ion temperature oscillates with the 
local plasma frequency. A phenomenological model was developed in [71] to 
describe the ion equilibration and kinetic energy oscillations in a manner that 
could extract information on damping. The main assumptions of this model 
are that ions execute damped harmonic motion at their local plasma oscillation 
frequency and Eq. (1631) determines the equilibrium ion kinetic energy and the 
amplitude of oscillation around this value. Fits to the data using this model 
showed that the damping time is close to l/^ P) i, and higher particle density and 
lower electron temperature give slightly stronger damping. Numerical studies 
of homogeneous strongly coupled plasmas [5Uf2"9"] suggest a similar damping 
time for Ti > 5, with lower T\ leading to faster damping. 



4-3 Collective Electronic Plasma Modes 

Oscillation of the electron plasma, or more precisely of the electron density, is 
a fundamental collective plasma mode first described by Tonks and Langmuir 
[138J. Excitation of this mode in an early experiment [32J firmly established 
that the photoionized ultracold gas forms a plasma. 

A straightforward way to resonantly excite this mode is with a radio-frequency 
(rf) electric field (Fig. [5]). A small DC bias directs electrons that escape the 
plasma to a charged particle detector in order to measure the response. Typical 
electron signals from such an experiment [32] are shown in Fig. [20h . The 
plasma is created at t — and some electrons promptly leave the sample and 
arrive at the detector at about 1 /j,s, producing the first peak in the signal. The 
resulting excess positive charge in the plasma creates the Coulomb potential 
well that traps the remaining electrons [2] . As the plasma expands, the depth 
of the Coulomb well decreases, allowing the remaining electrons to leave the 
trap. This produces the broad peak at ~ 25 /is. 
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Fig. 20. Electron signals from ultracold xenon plasmas created by photoionization 
at t = 0. (a) 3 x 10 4 atoms are photoionized and E e /kB = 540 K. Signals with 
and without rf field are shown. The rf field is applied continuously, (b) 8 x 10 4 
atoms are photoionized and E e /k-Q = 26 K. For each trace, the rf frequency in 
MHz is indicated, and the nonresonant response has been subtracted. The signals 
have been offset for clarity. The resonant response arrives later for lower frequency, 
reflecting expansion of the plasma; Reused with permission from |32j . Copyright 
2000, American Physical Society. 

In the presence of an rf field an additional peak appears in the electron signal, 
indicating that resonant excitation of plasma oscillations has pumped energy 
into the plasma and raised the electron temperature. The elevated electron 
temperature causes more electrons to reach the detector because it increases 
the evaporation rate of electrons out of the Coulomb well. 

This data [32] was initially analyzed assuming an expression for the angular 
oscillation frequency that is valid in a homogeneous gas (uj p>e = ^Je 2 p e /eom c 
[138J), and assuming that the excited mode was localized in regions of near res- 
onant density. By assuming the peak response arrived when the average den- 
sity in the plasma was in resonance, the changing resonant frequency showed 
the decrease in plasma density with time as the plasma expanded into the sur- 
rounding vacuum (Fig. [20b). Section I4~4l will discuss the expansion dynamics 
in detail. 

Reference jl9] pointed out that even in the zero-temperature limit, the elec- 
tron plasma resonance frequency should be modified because of the inhomo- 
geneous density distribution. This effect also causes the well known result 
that for a spherically symmetric flat-top density distribution, the frequency 
becomes lj = cu p , e / v3- For a Gaussian density distribution, the plasma has 
a continuous spectrum of electron-density oscillation frequencies, but spectral 
weight is concentrated in a damped quasimode that is concentrated in lower 
density regions of the cloud and has a frequency between 1/4 and 1/3 of u> p>e 
for the peak density in the distribution. This implies that the approach of [32] 
underestimated the density by about a factor of three, but does not change 
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Fig. 21. Electron emission from an expanding ultracold plasma. Data is taken as 
in Fig. [20] but the plasma is more dense and the applied rf field is more intense. 
For a given rf frequency, the earliest rf response peak is the cold plasma mode at 
the frequency predicted by the theory of [39]. The later modes are thought to be 
Tonks-Dattner resonances |139| . Reused with permission from [501. Copyright 2006, 
American Physical Society. 

the general conclusion that the plasma oscillation can be used to follow the 
evolution of the plasma density. 

In very recent work, this idea has been extended by Fletcher et al. [50J who 
observed what appear to be Tonks-Dattner modes, which are electron density- 
waves characteristic of plasmas with inhomogeneous density distributions in 
which thermal effects are important. (See Fig. [2TJ) If one goes beyond cold 
plasma theory, and includes the effect of electron temperature, electron plasma 
oscillations follow the Bohm-Gross dispersion relation |14Uj 



where k(r) is the local wavenumber. When the temperature or wavenumber 
becomes high, the thermal term becomes dominant and the electron density 
oscillation becomes a sound-wave. Tonks-Dattner modes are resonant electron- 
density standing waves resulting from this thermal effect, and they were orig- 
inally seen in cylindrical plasmas designed to simulate meteor tails |141lll42j . 
In these experiments, the plasma density increased towards the central axis, 




(64) 
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and the modes resonated between the outer edge of the plasma and the re- 
gion where u = w PjB . For higher frequency modes, more wavelengths fit in the 
resonant region. It is at first surprising that thermal effects are important in 
ultracold plasmas. But because of the small plasma size, the wavelength for 
density waves must be small and the wavevector must be large. This makes the 
thermal term in the Bohm-Gross dispersion relation (Eq. fl64|) ) comparable or 
larger than u Pie . 

An accurate quantitative description of Tonks-Dattner modes was given by 
[139J, and it was adapted to explain the observations in ultracold plasmas. 
The agreement is striking in spite of the fact that the authors had to make a 
somewhat ad hoc assumption about the outer boundary for the standing wave 
since there is not a well-defined edge to a Gaussian-shaped ultracold plasma. 
With additional theory input to understand this point, Tonks-Dattner modes 
should provide an accurate diagnostic of the electron temperature with good 
temporal resolution. 

Hence, both types of collective electron excitations can provide information on 
the time evolution of the plasma, its changing density and electron tempera- 
ture, and may be used to monitor dynamical effects in the plasma expansion. 



4-4 Plasma Expansion 



4-4-1 Physical Description 

In all experiments performed so far, the plasma is unconfined and expands into 
the surrounding vacuum. On a time scale of about 10 /xs, the plasma changes 
in size noticeably, as shown with absorption images of strontium plasmas in 
Fig. El 

As discussed in Section |4.1 .11 in conjunction with the initial electron trapping, 
the local charge imbalance due to the finite electron temperature causes a con- 
fining space charge potential. This traps the electrons such that a quasineutral 
plasma state develops, but the potential experienced by the ions has opposite 
sign, which drives the expansion of the plasma at later times. In an alter- 
native, but equivalent description, simple hydrodynamic arguments [32] show 
that the expansion can also be related to the thermal pressure of the electrons, 
since the hydrodynamic force per ion due to the electron thermal motion and 
plasma spatial gradient is given by (see eq. (l2T]) ) 



- Vpi = — — r . (65) 



Pi o" 2 



This yields the same results as found from collisionless plasma dynamics using 
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Fig. 22. False color images of expanding ultracold neutral strontium plasmas. The 
time after photoionization is indicated. 

the Vlasov equation as described in Section IXTTlR 

This model predicts that the plasma phase space density retains its Gaussian 
shape (Eq. fl22l) ) during the expansion, with a quadratically increasing rms-size 
(Eq.|MU) 

a(t) 2 = a(0) 2 + vlt 2 , (66) 



where the expansion velocity v = \J k^T e (fi) / rrii is solely determined by the 
initial electron temperature. The mechanical effect of the electrons on the ions 
is to induce a radially outward directed ion velocity (Eq. [26]) 



knTJO) t , , 

u r, t) = cl ' 2f2 r 67 

that is distinct from random thermal velocity, such that the ion kinetic energy 
due to expansion increases as E\ = \k-QT e {Q)v 2 t 2 / a 2 {t) . This energy increase 
has to be compensated by a decreasing temperature of the electrons which cool 
according to T e (t) = T e (0)a 2 (0)/a 2 (t) (Eq. MM- Tnus , asymptotically all the 
initial electron energy is transfered to the ions. Note that this energy transfer 
is not due to elastic collisions between the electrons and the ions, which would 
require milliseconds to equilibrate the energy of two plasma components. 

The expansion of neutral plasmas in different geometries has been intensively 
studied for decades, and much theoretical work has recently been done on the 



3 Strong coupling of the ions causes negligible modifications of the expansion 
dynamics (Section I3.2.3|) . 
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Fig. 23. Expansion of the plasma for N = 5 x 10 5 photoionized atoms. The expansion 
is described by n = N/[4ir(aQ + VQt 2 )] 3 ^ 2 , where o"o is the initial rms radius, Vq is 
the rms radial velocity at long times, and n is the density in resonance with the rf 
field, which is assumed to equal the average density in the plasma. Horizontal error 
bars arise from uncertainty in peak arrival times in data such as Fig. [20b. The fits 
were consistently poor at low E e , as in the 3.9 K data. Reused with permission from 
[32J. Copyright 2000, American Physical Society. 

dynamics of quasineutral plasmas expanding into a vacuum [9"2"|143|ll44lll45j . 
This problem is important for understanding plasmas created by irradiation 
of solid and thin film targets [146|17|16|19|20j . clusters [25], and gas jets [27] 
by intense laser pulses, and dusty plasmas |147j . which are important in as- 
trophysical systems and plasma processing. 



4-4-2 Studying the Expansion with Electron Plasma Oscillations 

The first quantitative experimental study of the expansion was performed 
using excitation of electron plasma oscillations as a probe of the density as 
described in Section l4~3l Data such as in Fig. [231 was fitted assuming a{t) 2 = 
cr(0) 2 + (vot) 2 for the expansion. The oscillation probe is only sensitive to 
times when the plasma has already reached a final velocity. Figure [23] shows 
the extracted values of v . For data with E e > 70 K, the expansion velocities 
approximately follow t> = [E e / (am;)] 1 / 2 , where m x is the ion mass and a = 1.7 
is a fit parameter. This equates to vq = [0.9fcsT e (0)/mi] 1//2 if fcsT e (0) = \E ei 
in reasonable agreement with the theory. 

For very low E e , the plasma expands as if the electrons had much more energy 
than E e (Fig. [21]). Different authors addressed this issue in a series of theory pa- 
pers pointing to the importance of continuum lowering [148] , disorder-induced 
heating [37f38] , and recombination |42j as sources of electron heating. A more 
detailed discussion of these issues is given in Section I4T51 The heating becomes 
important just as r e is approaching one, which is interesting, as it might in- 
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Fig. 24. Expansion velocities, vq, found from fits to data such as in Fig. [23j The 
initial average density, no, varies from 6 x 10 6 to 2.5 x 10 9 cm . The solid line, 
vq = \J E e / cmi\, with a = 1.7, is a fit to data with E e /ks > 70 K. The behavior 
of low E e data is discussed in the text. Uncertainty in vo is typically equal to the 
size of the symbols. There is a 0.5 K uncertainty in E e /kB reflecting uncertainty 
in the dye laser wavelength. Note that for E e /kB < 70 K, vq shows a systematic 
dependence on no- Reused with permission from [32]. Copyright 2000, American 
Physical Society. 

dicate the onset of strong coupling effects of the electrons. On the other hand 
this heating also implies that it is difficult to create strongly coupled electrons. 
There is still more work to be done to study the dependence of the heating on 
various parameters in order to work out the details of the dominant heating 
sources. 



4-4-3 Studying the Expansion with Optical Probes 

The expansion can also be studied with optical probes. Because of the better 
time resolution of these techniques, it is possible to follow the entire expansion 
from plasma formation, through ion acceleration, to terminal velocity. One can 
thus make an accurate comparison between experiment and theory. 

References [70] and [13] used the ion fluorescence to follow the spatial expan- 
sion of a calcium plasma. By placing the focused excitation laser at various 
distances from the center of the plasma, they measured fluorescence profiles up 
to 50 fis after photoionization (Fig. 1 of [13]) and found reasonable agreement 
with a modified version of the theory [33f73 J that addressed the lack of spher- 
ical symmetry. They also found that for E e < 100 K, the plasma expanded 
with more energy than expected from E e [70]. The ion velocity was probed 
spectrally using fluorescence from a resonant excitation beam passing through 
the plasma center [13]. Changes in the fluorescence level could be ascribed to 
changes in Doppler broadening and ion velocity, resulting in the data in Fig. 
M 
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Fig. 25. Studying plasma expansion with optical absorption imaging. The rms ve- 
locity along the laser beam is extracted from the width of the absorption image 
(|2.2.2p . The dominant contribution is from the expansion velocity but there is a 
contribution from the ion thermal velocity that produces the lOm/s offset at short 
times. The initial peak density of these plasmas was no = 3.5 x 10 15 m -3 and initial 
size <To = 1 mm. The initial electron kinetic energies are indicated in the legend. 
The solid-line fits to the data assume a self-similar expansion of the plasma as de- 
scribed in Eq. [67l This implies that the 1-dimensional rms ion velocity is given 
by y/kBT^ c fi/mi, where T^ e ff is given by Eq. (|B.4|) . The initial electron temper- 
ature can be extracted from the fits, and the agreement with expected values for 
the two higher temperature plasmas indicates that three-body recombination and 
disorder-induced heating of electrons is not significant for this data. If the electron 
temperature did not drop due to adiabatic expansion, the ion velocity would follow 
the dashed- line curves. 



In [28], the ion kinetic energy was found with the absorption spectra, and 
this was related to the expansion velocity and the average ion temperature. 
Initially, the technique was only applied to early times before adiabatic cool- 
ing of the electrons became important, but recent improvements have allowed 
study of longer times (Fig. [251) . For a broad range of conditions, agreement is 
excellent with Eq. flB.41) . which assumes a self-similar expansion (Eq. [21]) and 
that after the first few microseconds the only mechanism effecting the elec- 
tron temperature is adiabatic cooling. At lower E e and higher density, similar 
measurements show that electron heating effects contribute more significantly 
and for longer times (Section 14.5.21) . The data shows strong indications of 
the approach to a terminal velocity, connecting the acceleration phase to the 
measurements of Figs. [23] and [24] 
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Fig. 26. Spatial densities p\ (solid) and p a (dashed) of the ions and recombined 
atoms, respectively, at t = 3 ^s (a) and t = 31.3 us (b), compared to the Gaussian 
profile assumed for the hydro dynamical model (dotted). Additionally, pi obtained 
from the particle simulation using the mean-field interaction only is shown as the 
dot-dashed line in (a). Reused with permission from [73]. Copyright 2004, American 
Physical Society. 

4-4-4 Plasma expansion and ion dynamics 

The effect of the plasma expansion on the ionic temperature and, in turn, the 
influence of the ion relaxation on the expansion, can be understood qualita- 
tively from the hydrodynamical description of Section 13 . 11 First, neglecting 
ionic correlations (U# = 0) the expansion leads to adiabatic cooling of the 
ions, analogous to the electron dynamics (see Eq. ( 1241) ). In addition, however, 
the ion temperature also changes due to the development of ionic correla- 
tions mainly through disorder-induced heating discussed in Section 14.2.11 and 
described by eqs. (|39|) . Moreover, these correlations reduce the ion- ion inter- 
action (compared to the spatially uncorrelated case, where pairs of closely 
neighboring atoms carry a large potential energy), and therefore lead to an 
effective negative acceleration, expressed by the Ua/3-term in eq. (I39bl) . in 
addition to the ideal thermal pressure. This contribution, which corresponds 
to the average nonideal pressure known from homogeneous systems [31197] , also 
leads to an effective potential in which the ions move. As they expand in this 
potential, the thermal energy increases due to energy conservation. Finally, 
( I39bl) shows that the expansion is to the largest extent driven by the elec- 
tronic pressure, since T e ^> Tj even after the disorder-induced heating, since 
C/jj is of the same order as T ; . 

As the plasma expands, the spatial profile of the ions must deviate from its 
original Gaussian shape [33]. This is mainly due to deviations from quasineu- 
trality, e.g., deviations from the linear space dependence of the outward di- 
rected acceleration, at the plasma edge. The influence of the nonlinear corre- 
lation pressure on the density profile is of minor importance, as can be seen 
by comparing the solid and dot-dashed line of Fig. [26h in the inner plasma 
region. As known from earlier studies of expanding plasmas based on a mean- 
field treatment of the particle interactions [149,150 33], a sharp spike develops 
at the plasma edge, shown by the dot-dashed line in Fig. I26h . At later times, 
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Fig. 27. Measured time evolution of the electron temperature (a) and electronic 
Coulomb coupling parameter (b) for initial ion numbers of ~ 10 6 , initial plasma 
sizes of 0.4 mm< a < 0.7 mm and different initial electron temperatures. Reused 
with permission from [67] . Copyright 2004, American Physical Society. 



this spike decays again when the maximum of the hydrodynamic ion veloc- 
ity passes the position of the density peak, so that the region of the peak 
is depleted. Ultimately, at long times, the entire plasma approaches the pro- 
file of a quasineutral selfsimilar expansion |150j . From Fig. I26h it becomes 
apparent that with ionic correlations, the peak structure is less pronounced 
than in mean field approximation. This is due to dissipation caused by ion-ion 
collisions which are fully taken into account in the hybrid-MD simulation. As 
shown in [150j . by adding an ion viscosity term to the hydrodynamic equations 
of motion, dissipation tends to stabilize the ion density and prevents the occur- 
rence of wavebreaking which was found to be responsible for the diverging ion 
density at the plasma edge in the case of a dissipationless plasma expansion. 
Furthermore, the initial correlation heating of the ions largely increases the 
thermal ion velocities leading to a broadening of the peak structure compared 
to the zero-temperature case. 
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Fig. 28. Simulated time evolution of electron temperature scaled by the electron 
density with T e (0) = 0. Reused with permission from [41]. Copyright 2002, American 
Physical Society. 



4-5 Electron Heating Mechanisms 



Figures [23] and [24] suggest that plasmas strongly deviate from the the sim- 
ple selfsimilar expansion (Eq. (|24"1) ) when the plasma becomes dense and the 
initial kinetic energy of the electrons, E c , becomes small, suggesting that addi- 
tional physical processes heat the electrons. This question has been addressed 
in a series of theory papers [4~3"|148II38|I127|I41|I37II42 .33 . 1 5 1 J . proposing different 
heating mechanisms. The general conclusion is that the electron heating pro- 
cesses restrict the temperature of the electrons such that the electron Coulomb 
Coupling parameter is limited to r c <,0.2. 

The NIST group [HI] determined the electron temperature directly as a func- 
tion of time after photoionization by tipping the confining Coulomb potential 
and measuring the fraction of electrons that escaped over the lowered barrier. 
The potential was tipped by applying a small electric field pulse. The result 
in Fig. [27| shows that over a wide range of initial electron kinetic energies, 
E e , the electrons cooled or heated within a couple of microseconds to T e « 60 
K. T e subsequently dropped over the next 10 [is along a similar curve for all 
conditions. While the high-temperature data are consistent with the collision- 
less expansion (Eq. ([24]) ). the low-temperature results clearly indicate electron 
heating as the final temperature is well above the initial kinetic energy given 
to the electrons. 
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4-5.1 Disorder-induced electron heating 

While the temperature measurements of [67] were limited to times larger than 
a few microseconds due to the finite duration of the electric field pulse, detailed 
numerical studies of the short-time plasma dynamics have been performed by 
several groups [38|l41|37j . Since an accurate molecular dynamics treatment 
requires huge numerical efforts, in [38, 4Tf3~T] this aspect of the plasma dynam- 
ics has been studied for a smaller model system with a reduced ion-electron 
mass ratio. Fig. [28] shows the calculated time evolution of the scaled electron 
temperature 1/T e , for the extreme case of both zero initial ion temperature 
and zero initial electron kinetic energy, i.e., for an ionizing laser frequency 
that exactly matches the ionization potential of the cold atoms. Naively, one 
might expect very strong particle correlations at these low temperatures. How- 
ever, the electron component is not created in a correlated equilibrium state 
and heats up in the same way as dicussed in Section 14.2.11 for the ionic case. 
Immediately after photoionization spatial electron-electron and electron-ion 
correlations are practically negligible. They begin to develop on a timescale 
of the inverse electronic plasma frequency u~l = ^Jm e eo/e 2 p e lowering the 
electron potential energy, which is released as heat to the electrons. Due to 
the small electron mass the correlation build-up is much faster than for the 
ions, heating up the electrons within a few tens of nanoseconds. Since there are 
two contributions from electron-electron and from electron-ion correlations, we 
may estimate the final equilibrium excess correlation energy as two times the 
expression for u COTT given in Eq. (j4"3"]) . Energy conservation, 3T~ 1 /2 = 2w corr , 
yields T e pa 1, in agreement with the temperature dynamics shown in Fig. [281 
An estimate of the released energy for the expansion velocity measurement of 
Fig. [Ml yields a plateau velocity of 13m/s. The disorder induced heating alone 
is therefore not sufficient to explain the observed enhanced plasma expansion. 



4-5.2 Rydberg atoms formed through three-body recombination 

In three-body recombination (TBR) |152|T4"4"] , an electron and an ion recombine 
to form a highly excited Rydberg atom, while a second electron participates 
in the collision to conserve energy and momentum. TBR is an important pro- 
cess in ultracold neutral plasmas, which may be surprising because at typical 
plasma temperatures of 1000 K or more, TBR is only important in dense plas- 
mas. At low densities on the order of 10 10 cm -3 , recombination typically pro- 
ceeds by radiative recombination (RR) and dielectronic recombination (DR) 



TBR dominates in ultracold plasmas because the total three-body rate varies 
with temperature as T~ 9 / 2 . This divergent temperature dependence has mo- 
tivated investigations into whether TBR theory must be modified in the ul- 
tracold regime |154II155|I156||4"3"] . In addition to the fundamental interest in 



57 




50 100 150 200 

Time [/is] 



Fig. 29. Electron signal from a plasma created by photoionization of 10 5 atoms at 
t = 0, with E e /ks = 206 K. The first and second features represent free electrons es- 
caping from the plasma. The third feature arises from ionization of Rydberg atoms. 
A field of 5 mV/cm is present before the large field ramp commences at about 120 [is. 
Electrons from the first two features are collected and detected with 10 times less 
efficiency than the Rydberg electrons. Reused with permission from [66J. Copyright 
2001, American Physical Society. 

TBR at low temperatures, the process is important in the formation of cold 
antihydrogen through positron-antiproton recombination [4"5l4 6j. In ion stor- 
age rings, collisions between ions and electrons in electron coolers can be in 
the 10-100 K energy range, and anomalously large recombination rates have 
been observed for highly charged ions [157] . Some authors have suggested an 
enhanced contribution from TBR to explain a portion of the excess[158j. 

In addition to the fast initial rise of the temperature due to disorder-induced 
heating, Fig. [28] indicates further electron heating on a slower timescale, which 
is due to the formation of highly excited Rydberg atoms. Figure [29] shows 
experimental evidence of recombination in an ultracold neutral plasma. The 
origin of the time axis corresponds to the photoionization pulse. Immediately 
after that, some electrons leave the sample and are directed to the detector by a 
small DC field (5 mV/cm), producing the first peak in the signal. The ions are 
essentially immobile on this timescale, and the resulting excess positive charge 
in the plasma creates the Coulomb potential well that traps the remaining 
electrons [2]. The broad peak corresponds to the escape of free electrons as 
the Coulomb potential well vanishes as the plasma expands. After the free 
electrons have escaped, the electric field is increased to 120 V/cm in 50-100 fis. 
This field can ionize Rydberg atoms bound by as much as 70 K, corresponding 
to a principal quantum number of about n = 47. From the number of electrons 
reaching the detector, the number of Rydberg atoms formed can be inferred, 
and from the field at which the atoms ionize the distribution of Rydberg atoms 
as a function of n can be constructed. The Rydberg atoms survive for hundreds 
of microseconds with no significant change in their distribution. This implies 
that most of the atoms are in high angular momentum states which have long 
radiative lifetimes [58J. 



58 




Fig. 30. (left): Comparison of PICMCC-calculations [32] and measurements [32] of 
the plasma expansion velocity. Initial state parameters are the same as in fig|24i 
Reused with permission from [32]. Copyright 2002, American Physical Society, 
(right): Calculated time evolution of the electron temperature [29] compared to 
experiment [67] for two different initial electron temperatures of 6.7K (filled circles) 
and 66.7 (open circles). The remaining parameters are the same as in fig!27l Reused 
with permission from [29j. Copyright 2005, Institute of Physics. 

The classical three-body recombination theory, which is well confirmed for 
high temperature plasmas, predicts a total recombination rate of 

« 3.8 • 10- 9 T- 9 /V e 2 s- 1 , (68) 

where T e is given in K and the density p e is in cm -3 . For typical parameters 
of T = 9K and p e = 3 • 10 9 cm _1 this yields an extremely large recombination 
rate of K^f = 2.4 ■ 10 6 s" 1 . Hence, one would expect the plasma to almost 
completely convert into a Rydberg gas within a few microseconds. Contrary 
to this, only a few percent of recombination is observed after 12/is for the 
above initial conditions [66]. This has led to a series of theoretical investiga- 
tions questioning the validity of the standard recombination theory [43] and 
proposing strong coupling effects on the recombination rate [154jll59| . Sub- 
sequent theory and experiments |66)37II42|[67] indicate the opposite, namely 
that the plasma is driven out of the strongly coupled regime, by subsequent 
electron-atom collisions involving the formed Rydberg atoms, which tend to 
additionally heat the electron gas. The corresponding increase in temperature 
is in fact more rapid than from direct recombination and quickly quenches fur- 
ther recombination due to the strong T e ~ 9//2 dependence of the rate. The simple 
hydrodynamic description with standard collision rates [44] . as described in 
Section I3.1.2[ yields a surprisingly good description of the plasma dynamics. 
As shown in Fig. [30] the enhanced expansion velocity of Fig. [23] can indeed be 
reproduced upon inclusion of electron-Rydberg-atom collisions, and the hy- 
drodynamic model even yields a quantitative description of the temperature 
measurement [67] . 

The measured distribution of principal quantum numbers for the Rydberg 
atoms also revealed that moderately tightly bound levels (n < 50) were popu- 
lated in the plasma, and that the binding energy of the Rydberg atoms closely 
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matched the excess energy that appeared in the expansion as shown in Fig. 
[24] [66]. 

Although the discussion about the validity of the conventional recombination 
theory under the conditions of an ultracold plasma is still not settled, the pre- 
sented measurements [66.67J together with supporting theoretical calculations 
provide strong arguments for the applicability of the classical recombination 
rates even in the low -temperature regime. 

4-6 Coulomb Coupling Parameters 

As stated many times before in this review, fundamental interest in ultracold 
neutral plasmas to a large extent stems from the possibility of creating strongly 
coupled plasmas. For a system of charges embedded in a uniform neutralizing 
background, formally called a one-component plasma (OCP) [HZ], local spa- 
tial correlations characteristic of a strongly coupled fluid appear for T > 2. As 
r further increases and eventually exceeds a value of 174 the system crystal- 
lizes. While the OCP model was originally developed to describe dense neutron 
stars, with ions embedded in a degenerate electron background, a clean exper- 
imental realization of an OCP is provided by laser cooled ion clouds confined 
in Penning [16U|ll61|l97j or Paul [162|I163] traps. 

In a real two-component plasma (TCP) the situation can be very different. 
Besides the possibility of forming bound states, electrons will tend to shield the 
ion-ion interaction potential. For weakly coupled electrons, within the Debye- 
Hiickel theory, the Coulomb interaction is replaced by an interaction potential 
V\\ = exp(— t/Ad) of Yukawa-type. This shielding can drastically reduce the 
ion-ion interaction energy and the correlations, as approximately expressed by 
an effective coupling constant T* = T^ K [3711134] . where k = a/X D = y/STl. 
While this is certainly a good approximation for T e -C 1 it breaks down as 
k > 1. 

Ultracold plasmas provide the first precise experimental studies of two-component 
plasmas with both components in or near the strongly coupled regime. As we 
have discussed, however, various heating mechanisms currently preclude the 
exploration of plasmas deep in the strongly coupled regime (r 3> 1). 

4-6.1 Electronic Coulomb Coupling Parameter 

At the lowest possible initial electron temperatures of about O.IK, limited by 
the linewidth of the ionizing laser, and typical densities of 10 9 cm -3 very large 
Coulomb parameters of about 30 might be expected. As outlined in the pre- 
vious section there are however various heating mechanisms quickly decreas- 
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Fig. 31. Time evolution of the electron Coulomb coupling parameter for an initial 
average density of 10 15 cm -3 and two different initial electron temperatures of 10 K 
(solid lines) and 50 K (dashed lines). Two simulations, including and excluding 
Rydberg atom formation are shown for each parameter set. Reused with permission 
from [42J. Copyright 2002, American Physical Society. 

ing this value. After an initial disorder induced heating phase, the electrons 
are further heated by Rydberg atom formation and subsequent electron-atom 
collisions. Interestingly, T e approaches an almost universal value largely in- 
dependent of the initial plasma parameters. This is demonstrated in Fig. [3T] 
showing results of PICMCC simulations (see Section 13. ip for two different 
initial temperatures [42]. When three-body recombination is turned off, T e 
increases in time due to the adiabatic cooling of the electron gas. When the 
formation of Rydberg atoms is taken into account this increase is strongly 
suppressed by the heating of the electron gas. While the initial evolution is 
qualitatively different for the high and low temperature case, respectively, 
both curves approach a value of r e « 0.2 within about 10/xs. This prediction 
has been confirmed by subsequent temperature measurements, discussed in 
Section 14.51 Although such a PICMCC treatment can not describe any type 
of particle correlations of the plasma and hence neglect the disorder-induced 
electron heating, the measured value of the plateau is consistent with the the- 
oretical result of 0.2, within the experimental uncertainty. This suggests that 
Rydberg atom formation indeed constitutes the dominant heating mechanism 
and correlation effects do not play a significant role for the electron dynamics. 

4-6.2 Ionic Coulomb Coupling Parameter 

This situation is different for the ionic component, whose temperature is prac- 
tically not affected by inelastic collisions due to the much larger mass of the 
ions. Heating can therefore only arise from correlation build-up. PICMCC 
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simulations (Fig. |3~T1) are unable to describe Tj, since an accurate treatment 
of ion-ion correlations is required, as can be achieved within the hybrid-MD 
approach introduced in Section [3.2.31 

In order to calculate T; one needs to know the ion temperature and den- 
sity. Experimentally, T; iave can be determined from the absorption spectra 
[28] or ion fluorescence [13J. The density can be determined from absorption 
images, knowledge of the neutral atom parameters and the ionization frac- 
tion, or the timescale of the disorder- induced heating. In order to estimate 
the effect of screening, the electron temperature must be known. Finally, 
the inhomogeneity of the density distribution must be taken into account. 
One way to do this is to calculate the density-averaged coupling parameter, 
r* javg = (exp[— K(r)]e 2 /[47T£ a(r)/c£TJ jave ]). Measurements find for strontium 
plasmas r* ~ 2 after the disorder-induced heating phase for a wide range 
of initial conditions [2B]- With lower T e and higher density, r* ave was only 
slightly higher. A similar result, r ~ 4, was found for calcium plasmas [13J. 

The surprisingly small variation in r* suggests that disorder-induced heat- 
ing is a natural feedback mechanism that leads to equilibration into a state just 
barely in the strongly coupled fluid regime independent of the initial condi- 
tions. This conclusion is supported by the qualitative discussion on disorder- 
induced heating in Section 14.2.1] where it was argued that the final tem- 
perature should be of the order of the Coulomb interaction energy between 
neighboring ions, which directly translates into r ~ 1. 

The question remains of how the Coulomb coupling parameter evolves in time 
as the system expands. The simple hydrodynamical model (Eq. [241) predicts 
that the ion temperature decrease during the expansion should overcompen- 
sate the accompanying decrease in density, so that the Coulomb coupling 
parameter should increase oc p -1 / 3 without bounds as the system expands. 
Experimentally, such a prediction is hard to check, since the directed radial 
expansion velocity quickly overwhelms the thermal velocity component. The- 
oretically, on the other hand, it is straightforward to separate the thermal 
velocity components in the hybrid-MD simulation |39j . 

Figure [32] shows that 1^ does increase as the system expands. Hence, one 
might conclude that it is sufficient to simply wait long enough to allow the 
system to cool down in order to reach the crystallized state. However, it is 
questionable whether one should speak of a "crystal" if the density drops 
below that of interstellar space by the time it reaches a T of 200. But there is 
also a fundamental problem; the coupling parameter measures the amount 
of correlations in a plasma if the system is in in thermodynamic equilibrium, 
but the plasma is created far from equilibrium in a completely uncorrelated 
initial state. Disorder-induced heating leads to local local thermal equilibrium, 
but the plasma continuously expands, so a true static equilibrium is never 
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Fig. 32. Ionic Coulomb coupling parameter for a plasma with iVi(O) = 5 • 10 4 , 
Pi(0) = 1.1 • 10 9 cm~ 3 and T e (0) = 50K. The solid line shows the Coulomb cou- 
pling parameter calculated from the average temperature and density, the dashed 
line marks the coupling parameter extracted from pair correlation functions (see 
text). Inset: blow-up of the short-time behavior. Reused with permission from [30] . 
Copyright 2005, American Physical Society. 

reached. One must compare the coupling parameter formally calculated from 
temperature and density with the spatial correlation properties of the system 
at all times. 



In [30], this question was studied by extracting an alternative coupling param- 
eter Ti from the pair correlation functions gn(r) obtained with the hybrid-MD 
approach (dashed line in Fig. 1321) . In line with the preceding discussion, Ti 
and Ti initially have no similarity at all, reflecting the fact that the plasma is 
very far away from any kind of equilibrium. (In addition, Fig. 3(a) of reference 
[30] shows that the pair correlation function at early times is not well fit by 
that of an equilibrium system, so that a parametrization of by a single 
scalar quantity F\ does not adequately describe the plasma's spatial correla- 
tion properties. For 5 < u> Pt it < 20, the two curves for T; and T; approach each 
other and run in parallel, suggesting that the plasma is close to a changing 
quasi-steady state during this stage. 

For u Pti t > 20, however, the two curves diverge again, and f\ settles towards a 
finite value of Ii « 11 deep in the strongly coupled regime, but still far away 
from a value where the occurrence of long-range order might be observed. 
In other words, while the system expands and cools down (reflected in the 
further increase of 1^), the spatial correlations in the system (parametrized by 
r ; ) do not change anymore. On a formal level, this freeze-out of correlations 
is caused by another violation of the Bogoliubov assumption [136] . namely a 
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cross-over of the timescale for relaxation of correlations with the hydrodynamic 
timescale on which macroscopic system parameters change. More intuitively, 
the build-up of spatial correlations requires strong interactions between the 
ions in order to cause a rearrangement of their spatial correlational structure. 
As the system expands, the timescale on which such rearrangement takes place 
becomes longer and longer as the inter-ionic forces become weaker and weaker. 
At some point, the system simply becomes so dilute that the rearrangement 
does not occur anymore, and the spatial structure of the plasma freezes out 
in a selfsimilar expansion. As a consequence, the plasma remains in a state 
governed by strong short-range correlations, but no long-range order should 
form. 



5 Achieving Strong Coupling 



As emphasized in various sections of this article, one of the main motivations 
for studying ultracold plasmas from a plasma-physics perspective is the fact 
that they are strongly coupled, and one might hope to create crystalline plas- 
mas. Coulomb crystallization is actively studied in one-component systems, 
such as magnetically confined ions [16411165] and dusty plasmas [16 6|167lll68j . 
Equivalent studies in two-component plasmas, on the other hand, are more 
scarce, as their experimental realization poses several additional constraints 
[169J. In this respect ultracold neutral plasmas have attracted considerable 



attention as a promising opportunity to observe this phenomenon in the lab- 
oratory. 

However, as discussed above, three-body recombination strongly heats the 
plasma electrons, while the disorder-induced heating rooted in the (spatially) 
uncorrelated nature of the initial plasma state (i.e., after creation of the plasma 
by photoionization) rapidly heats electrons and ions. This was realized soon 
after the first experiment, and different suggestions have been put forward to 
mitigate these effects. 

For ions, these suggestions may be divided into two groups. On the one hand, 
one may avoid the disorder-induced heating. Since this heating effect is directly 
related to the uncorrelated initial state of the plasma, such approaches always 
aim at introducing some amount of correlation in the spatial distribution of 
the atoms before the plasma is formed by photoionization. We will discuss such 
approaches in Section 15.1.21 below. On the other hand, one might counteract 
the disorder-induced heating by removing energy from the system during its 
evolution rather than avoiding the heating from the outset. In this case, one 
needs a way to couple the system to its environment in order to dissipate part 
of its internal energy, i.e., to cool it. Doppler laser cooling has been suggested 
as a way to achieve this goal. Since it is, in our opinion, the most promising 
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perspective for achieving plasma crystallization at this point, we will discuss 
it in detail in Section 15.1.11 Finally, a recent proposal to increase the coupling 
of the plasma electrons by adding Rydberg atoms to the plasma, the only 
suggestion so far aiming at the electronic component of the plasma, will be 
briefly discussed in Section 15.21 

5.1 Ionic coupling 

5.1.1 Laser-cooling the ions 

As stated above, the idea behind laser cooling the plasma ions is to coun- 
teract disorder-induced heating by removing energy from the system via an 
external cooling mechanism. Such an approach is experimentally feasible if 
alkaline earth elements are used for the plasma, as their singly charged ions 
have a single outer electron which provides a convenient cooling transition for 
currently available laser systems. Hence, a "standard" Doppler cooling scheme 
may be used to cool the ions, with an optical molasses created by six coun- 
terpropagating lasers in a 3d configuration (cf. Fig. [2]). The plasma ions can 
then be cooled during the whole evolution of the system. 

Experimentally, the scheme outlined above has not yet been realized. However, 
there is little doubt that it should be feasibl^]. In fact, a first step in this 
direction has already been taken with the Doppler imaging of the ion dynamics 
described in detail in Section 12.2.21 since the same transition that is used for 
the imaging can also be used for the cooling. Theoretically, the above scenario 
has been studied in detail in [S^IIMj . 

Originally, the idea of Doppler laser cooling the ions was proposed in [37] 
and |17U] . These studies also correctly pointed out the potentially crucial role 
of elastic electron-ion collisions in the plasma in such a scenario. In a freely 
evolving plasma, such collisions are usually negligible, due to the fact that the 
timescale for electron-ion thermalization is much longer than the timescale 
of the plasma expansion. Hence, the heating of the ions during the first mi- 
crosecond of evolution due to these collisions is limited to a few milli-Kelvin, 
which is negligible compared to the additional disorder-induced heating which 
is orders of magnitude stronger. In the case of laser cooling, however, the final 
ion temperature one is aiming for is also in the milli-Kelvin range, i.e., of the 
same order as the energy transferred from electrons to the ions via elastic colli- 
sions. Hence, these collisions should be taken into account when discussing the 
possibility of plasma crystallization. In [3T|17U] , their influence on the achiev- 
able ionic temperature was estimated using reasonable values for initial-state 

4 It should be noted, though, that the required cooling rates are more easily 
achieved for light ions such as Be. 
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temperatures and densities. It was concluded that the ionic Coulomb coupling 
parameter would be limited to a value T; m 50 due to these elastic collisions, 
outside of the crystallized regime. However, these estimates assumed a con- 
stant heating rate corresponding to the chosen initial conditions. What was 
not taken into account is the expansion of the plasma, which strongly reduces 
the electron temperature and density, which reduces the collisional heating. 
As a consequence, these estimates are too pessimistic. 

A detailed numerical investigation of the laser cooling process was carried 
out in |34,94|, using both the hydrodynamic approach and the hybrid-MD 
approach introduced in Section [31 Insight into the modifications induced by 
the laser cooling can be obtained from a hydrodynamical description. On the 
level of a kinetic description, the influence of the cooling laser can be modeled 
via a Fokker-Planck term in the time evolution of the ionic distribution fr, 




V v (v/0 + ^A v/i 



(69) 



where (3 is the cooling rate and Tp the Doppler temperature. A derivation 
analogous to Section 13.11 leads to the macroscopic equations for the evolution, 



J^ 2 = 2 7( x 2 (70a) 



d_ _ (k B T c + fc B Ti 
dt m\a 



1= ^ B l,r i; -7(7-/3) (70b) 



JkfcuTO = -2 7 M - 2(3 (T, - T c ) (70c) 
at 

^{k B T c ) = -2 1 k B T e . (70d) 

As readily observed from Eq. (1701) the laser cooling dramatically modifies the 
expansion, not only quantitatively, but even qualitatively. While the width a 
of the plasma cloud increases linearly in time for the free expansion, it only 
varies as oc t 1//4 in the laser-cooled case. As a consequence, the plasma ex- 
pands much slowly than in the case of free expansion. This raises the hope of 
achieving a long-range ordered state, since we have seen in the previous section 
that it is the decrease in density and the corresponding increase in relaxation 
timescale that leads to a freeze-out of spatial correlation properties for a freely 
expanding plasma. In the laser-cooling scenario, the much slower decrease in 
density should allow for much longer times during which efficient rearrange- 
ment processes may occur, a prerequisite for the build-up of long-range order. 
One may, at first glance, think that this problem might be circumvented in 
any case by just increasing the cooling rate (3 sufficiently. However, this is not 
the case, since a too rapid cooling could freeze the plasma in a metastable 
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Fig. 33. Left: Distribution of scaled ion distances after an expansion time satisfying 
u!p t it = 240 (dots). The comparison with the correlation function of an OCP at 
T; = 700 indicates the development of strong order. The inset shows the evolution 
of the ionic Coulomb coupling parameter determined from the correlation energy 
(solid line) and from the ion temperature and density (dashed line). Right: Shell 
structure formation during the plasma expansion demonstrated by the radial ion 
density and a two dimensional cut through the plasma. Reused with permission 
from [34J. Copyright 2004, American Physical Society. 

disordered state. The system needs a certain time for spatial rearrangement 
to occur, which will not happen if the cooling is too fast. 

In order to obtain information about the build-up of spatial correlations, one 
needs to go beyond a kinetic description of the ion distribution function. In 
[M] this was achieved within the hybrid-MD approach of Section I3.2.3[ where 
the influence of the cooling lasers on the ions is described by a Langevin 
term consisting of a dissipative cooling force and a stochastic diffusive part. 
In addition, elastic electron-ion collisions should be taken into account in the 
description as discussed above. Such collisions can be implemented using a 
Monte Carlo procedure described in detail in [91]. It is found that the resulting 
structures formed during the system evolution may depend strongly on the 
initial-state parameters, such as number of ions, electron temperature, and 
cooling rate. Thus, the electronic temperature must be chosen sufficiently high 
in order to keep three-body recombination at a negligible level. As discussed in 
[M] , three-body recombination creates doubly-excited and hence autoionizing 
states with an excess energy of several thousand Kelvin above the continuum 
limit. Hence, any autoionization process results in a fast electron which would 
not be bound in the ionic potential well, thereby destroying the plasma. 

There are also significant constraints on the cooling rate f3. As discussed above, 
if the cooling is too fast, the plasma will be frozen in a disordered glass-like 
state. On the other hand, if the cooling is too slow, the density decreases too 
strongly for effective inter-ionic interactions to take place. As we have seen 
in Section I4.6[ reaching low temperatures and correspondingly high coupling 
parameters does not yet ensure spatial order in the plasma. However, it was 
demonstrated in [34"f94"] that, depending on particle number and cooling rate, 
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Fig. 34. Ionic intra-shell structure of the 1st, 3rd and 5th shell (from left to right). 
Reused with permission from [29], see also |34| . 

glass-like order or the formation of concentric shells may occur (see Fig. [33] 
and FigED. 

While the precise dynamics of this structure formation is not yet completely 
understood, it appears sufficiently robust that the crystallization should be 
observable experimentally with current state-of-the-art cooling techniques. 
Such an experimental confirmation, together with further detailed experimen- 
tal and theoretical studies, would allow significant physical insight into the 
dynamics of strongly coupled two-component plasmas at the liquid-crystal 
phase transition. In particular, there are both similarities and differences with 
laser-cooled ion clouds in ion traps, where shell structures are also observed 
and intensively studied [160f 171|I172II173|I174|I 17511176] . It may appear surpris- 
ing that long-range order may occur in the present case of an unconfined 
plasma where, in contrast to trapped ion clouds, the crystallization process 
can only take place over a finite amount of time. Pronounced differences also 
exist in the dynamics of the crystallization process, which starts at the outer 
edge of the ion cloud in the case of ions in a trap, but is observed to proceed 
from the center of the plasma in the present scenario. Moreover, the spherical 
structure of small ion crystals can be attributed to the spherical symmetry 
of the trapping potential, while in the present case it must be traced to the 
spherical symmetry of the initial density distribution. This raises the question 
of what is the expected symmetry resulting from an ellipsoidal initial-state 
density [13]. 

5.1.2 Other Approaches 

The crystallized regime would be reached more directly by avoiding the disorder- 
induced heating right from the start, rather than counteracting it by external 
cooling. Avoiding (or at least weakening) it requires a certain amount of spa- 
tial order in the initial state, i.e., one must generate a significant degree of 
spatial correlation before the plasma is created by photoionization. A number 
of proposals exist for how such an ordered state may be obtained, namely (i) 
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using fermionic atoms cooled below the Fermi temperature, so that the Fermi 
hole around each atom prevents the occurrence of small interatomic distances 
[10]; (ii) an intermediate step of exciting atoms into high Rydberg states, so 
that the interatomic spacing is at least twice the radius of the correspond- 
ing Rydberg state [103J; and (Hi) using optical lattices to arrange the atoms 
[103|Tl 33j . Calculations show that, indeed, these proposals could achieve a sig- 
nificant reduction of the disorder- induced heating, and the latter may even 
result in correlation-induced cooling of the plasma [133J. However, realizing 
the crystallized regime in this way turns to be very challenging experimentally 
and the laser-cooling approach discussed in the previous section seems to be 
more practical with current experimental capabilities. 



5.2 Electronic coupling 

All the scenarios discussed in the previous Section 15.11 aim at increasing the 
coupling of the ions in order to reach the crystallized regime. On the other 
hand, it would also be highly desirable to achieve strong coupling of the plasma 
electrons. Many previously employed assumptions, such as the dependence of 
three-body recombination rates on density and temperature, or the Thom- 
son cut-off for the existence of Rydberg states are expected to fail in this 
regime. Moreover, this would provide an intriguing realization of a classical 
two-component system where both components are strongly coupled, with 
exciting prospects for observing new and interesting phenomena. 

The only proposal so far targeting the electronic coupling parameter was in- 
troduced in [68J . The basic idea consists of deliberately adding Rydberg atoms 
to the plasma, which should then be collisionally ionized by the free plasma 
electrons, thereby reducing the kinetic energy of the electrons and cooling 
them. Such a scenario has been studied systematically in |177j . and it was 
found that a moderate but still significant manipulation of the electron tem- 
perature is possible in this way. Figure [35] shows the time evolution of T e for 
a plasma when Rydberg atoms with varying principal quantum number n are 
added at t = 0. As can be seen, adding Rydberg atoms with n = 36 to the 
plasma leads to a significant reduction of T e compared to the pure plasma 
without additional Rydberg atoms, in agreement with the original proposal 
[68] . On the other hand, this effect strongly depends on the principal quan- 
tum number of the added Rydberg atoms. For lower excitation (n = 26 in the 
present case) , even the reverse effect of a heating of the plasma electrons rather 
than a cooling can be observed. This behavior, which might be surprising at 
first glance, is connected with the existence of a bottleneck in the transition 
rates for the inelastic electron-impact induced bound-bound and bound-free 
transitions of the atoms. If the binding energy of an atom is larger than 
a critical (temperature-dependent) energy E c « 4k B T e , the atom is predom- 
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Fig. 35. Time evolution of T e for a plasma of Ni = 80, 000 cesium ions with 
T e (0) = 50 K and cr(0) = 155 fjm, after addition of iV R = 80,000 Rydberg atoms 
with principal quantum numbers n = 26,31,36 (dashed lines), compared to that 
of the pure plasma without additional Rydberg atoms (solid line). Reused with 
permission from [68j. Copyright 2005, American Physical Society. 

inantly driven to even lower states by collisions with the plasma electrons, 
thereby heating them. If, on the other hand, is smaller than E c , collisions 
tend to excite the atom to higher bound states or ionize it, and the plasma 
electrons are cooled as originally envisioned. 

Hence, by varying the principal quantum number of the added Rydberg atoms, 
it is possible to control the rate and even the sign of the change of the electron 
temperature. On the other hand, varying their number influences the absolute 
magnitude of the temperature change as well as the timescale over which 
cooling can be achieved. Together, a remarkable amount of control over the 
time evolution of the electronic temperature appears possible, in particular 
if one considers the possibility of more sophisticated schemes of adding the 
Rydberg atoms. Rather than adding all atoms simultaneously at t = 0, one 
may think of adding them continuously over a certain amount of time, or in 
a number of individual pulses (or generally speaking "on demand"). One can 
also change the principal quantum number of the atoms in time. (The latter 
might be necessary, e.g., to stay away from the bottleneck separating cooling 
from heating, which changes as T e changes in time.) 

In addition to the temperature reduction discussed above, the collisional ion- 
ization of Rydberg atoms also increases the density of the plasma. Thus, both 
effects act cooperatively, increasing the Coulomb coupling parameter T e as 
desired. It was found in |177j that, using relatively simple schemes for the 
addition of the atoms, T e reaches values of T e m 0.5, as compared to T e < 0.2 
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found in the experiments with pure plasmagjj. In more advanced scenarios 
as outlined above, reaching the strongly coupled state of T e > 1 should most 
probably be realizable. The coupling strengths obtainable are certainly far 
away from the regime where a crystallization of the electrons might be ob- 
served. Still, pronounced liquid-like correlation effects may be anticipated, 
and the first realization of a two-component system where both the electronic 
and the ionic component are strongly coupled promises exciting new develop- 
ments. 



6 Conclusions and Future Directions 

In this article we have reviewed the development in research on ultracold neu- 
tral plasmas since their first creation in 1999. We have described numerous ex- 
perimental tools to probe various aspects of the system behavior, ranging from 
optical imaging of individual ions to employing collective plasma excitations 
as probes for macroscopic quantities. By combining such techniques, tradition- 
ally rooted in both atomic and plasma physics, these experiments have studied 
a number of plasma physics problems. They have revealed several surprises, 
which fosters continuing theoretical interest. We have outlined the theoretical 
concepts, developed to describe and understand the non-equilibrium dynamics 
of ultracold neutral plasmas. Our discussion has covered diverse approaches, 
such as simple reduced dynamical equations, kinetic or hydrodynamic descrip- 
tions and demanding microscopic approaches. 

Based on the techniques introduced we have given a detailed description of the 
underlying physical processes driving this quite complex dynamical behavior. 
Within a hand-in-hand experimental and theoretical discussion a step-wise de- 
veloped physical picture of the plasma dynamics over various evolution stages 
from its creation until its complete expansion into the surrounding vacuum has 
been given. Thereby we have covered a wide spectrum of phenomena, ranging 
from low-temperature collisions and Rydberg atom formation to the physics 
of collective plasma excitations and effects of strong correlations between the 
plasma constituents. 

This wealth of interesting physics revealed in the past years demonstrates the 
rapid advancement in our understanding of these fascinating systems, but also 
raises various open questions, which we have pointed out in this article. Future 
studies in this young research field will certainly produce unexpected results, 
but several trends are likely to continue for the field. 



Experimental results so far are consistent with this theoretical prediction, but 
still inconclusive since current schemes of extracting an electron temperature from 
the experiment are not very precise. 
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As discussed in Section |5j the pursuit of stronger coupling in both the electrons 
and ions will remain an exciting and therefore intensively followed avenue of 
research. Experimental routes for increasing the ion coupling constant, such 
as laser-cooling [5711170] or photionizing a pre-ordered sample [1011133] . may 
enable the formation of Wigner Crystals [34J. This will hopefully allow ex- 
ploration of the phase diagram of particles interacting through a screened 
Yukawa potential [178|ll01|ll02j that will compliment studies of dusty plas- 
mas Creation of colder and denser electrons will show whether strongly 
coupled electrons can be realized in ultracold plasmas or three-body recombi- 
nation precludes this possibility 

Exploring the regime of stronger electron coupling will necessarily probe the 
complicated evolution of the electron temperature when rapid thermalization, 
disorder-induced heating, threshold lowering, and three-body recombination 
play significant roles. Powerful theoretical tools have been developed to address 
these phenomena [4Tf3~7 38,42,33,43,73]. but there is still need for a compre- 
hensive study of evolution in various regimes of initial conditions and a disen- 
tangling of the many contributing effects. As emphasized in [33], the synthesis 
of ideas from atomic physics and plasma physics is crucial for this research. 
While techniques exist to study electron temperature at later times [pTfoD] . 
at present there is no experimental diagnostic to cleanly explore the electron 
equilibration phase because electrons evolve so quickly. The development of 
such diagnostic would be a significant advance. 

As the field grows, there is great possibility to increase the connection to tra- 
ditional areas of plasma physics. The role of strong-coupling and the equilibra- 
tion and expansion dynamics partly resemble the behavior of laser-produced 
plasmas in high energy-density experiments. But as we have emphasized from 
the start of this review, the slower time scales, excellent diagnostics, and con- 
trol over initial conditions in ultracold plasmas provide unique opportunities 
to make significant contributions to these other areas. 

The recent success to produce cold Rydberg gases and plasmas in strong mag- 
netic fields of several Tesla [179|I180] has opened up a new regime, where the 
atomic scale and macroscopic system dynamics is fundamentally different from 
the field-free situation. Theory has just started to treat collisional [181II1U6|I107] 
and radiative |182|I183II184] processes under such conditions, while the macro- 
scopic dynamics of magnetized, ultracold, /xm-scale neutral plasmas still awaits 
a deeper understanding - all being essential questions for ongoing efforts to 
produce antihydrogen atoms [4oT4B] . 

New or improved experimental capabilities wait just over the horizon and 
will certainly open many avenues of research. As the resolution for imaging 
techniques approaches ~ 10 /im, this will allow studies of collective modes 
or other features of the ion density distribution. The use of light scattering 
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to directly observe spatial correlations [161] and the application of magnetic 
fields to confine the plasma will also provide new possibilities. 



In light of the rapid development reviewed in this article, the study of ultracold 
neutral plasmas promises to remain an exciting and growing field of research 
for many years to come. 



A Important Quantities 

CCP - Coulomb Coupling Parameter 
LTE - Local Thermal Equilibrium 
MD - Molecular Dynamics 
OCP - One Component Plasma 
TCP - Two Component Plasma 
PIC - Particle in Cell codes 

PICMCC - Particle in Cell Monte Carlo Calculation 
TBR - Three body recombination 
RR - Radiative recombination 
DR - Dielectronic recombination 

* in all expressions a stands either for ions (a = i) or electrons (a = e) 

particle number N a 

plasma temperature T a 

plasma width a 

plasma density p a 

peak density p = N a /(2ir<j 2 ) 3 / 2 

average density p — N a / (47rx 2 ) 3//2 

plasma frequency u Pja = [e 2 p a /(m Q ,e )] 1/ ' 2 

Wigner-Seitz radius a a = [3/(47rp a )] ly/3 

Coulomb coupling parameter T Q = e 2 /(47re a Q ,/i; B T), 

Debye screening length [£okBT a /(e 2 p a )] 1 ^ 2 

Bruckner parameter r s = a a / a® 

(ratio of interparticle distance to Bohr radius) 



B Extracting Plasma Parameters from the Absorption Spectrum 

The Voigt profile for the ion absorption cross section reads 



0, 7]) = J ds 



3*A 2 70 1 1 

\-y aH /2TTj 



(S - V Q ) 



2" 



a 



exp 



2a D (T i y 



(B.l) 
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where the parameters are described below Eq. (TT2T) in Section f2. 2. 21 



The expansion of the plasma described in Section 14.41 eventually becomes the 
dominant source of Doppler broadening for the ion absorption spectrum in an 
ultracold neutral plasma (Fig. [7]). Fortunately, the expansion is approximated 
well by a self-similar expansion of a Gaussian density distribution, and the 
dynamics takes a particularly simple form when there are no heat sources for 
the electrons and ion-electron thermalization is negligible. (See Section [3.1.21 
and [32.33.34j.) These assumptions are reasonable for plasmas of low density 
and high initial electron energy E e , although the exact regime of validity is 
the subject of current study. In this regime, the dynamics is determined by the 
parameters T e (0) and a, where T e (0) is the electron temperature after any ini- 
tial heating (Section 14. 51) . To be more specific, we assume that electron heating 
is only important at very early times before plasma expansion is important, 
and subsequently the electron temperature evolves only under the influence of 
the adiabatic expansion. Hence, /c B T e (0) may be slightly higher than \E e . 

The expansion velocity profile u(r, t) (Eq. (1261) ) gives rise to an average Doppler 
shift of the resonant frequency that varies with position, 



r / \ n rk B T e (0)t 

ov{r)=u z \ = — — cos# (t>.2) 

where r cos 9 is the displacement from the center of the cloud along the direc- 
tion of laser propagation and t is the time after photoionization. The position- 
dependent Doppler shift must be included in the exponent of the Gaussian 
describing the Doppler broadening in the Voigt convolution of the absorp- 
tion cross section a(i>, TJ) (Eq. fIB.lj) ). The latter becomes explicitly position 
dependent and takes the form 



a = Fir 

2vr 7 efr V 1 

F(r) = f ds- 



(B.3) 



l + 4f^_ N ) V27ra D (Ti(r)) 



exp 



[s -u - 5u(r) f 
2a D (T i (r)) 2 



where we have also included position dependence if the local thermal equi- 
librium temperature (Eq. fl63l) ). This variation of the ion temperature with 
position is the major complication to calculation of the full absorption S(u) = 
J <i 3 rpi(r)a(i/, 7], r). But by replacing T;(r) with the average ion temperature 
in the plasma, T ijave = / d 3 rTi(r) [54], and utilizing the spherical symmetry 
of the plasma, the integral over space can be evaluated, and the ion temper- 
ature can be replaced with an effective ion temperature (T i e fj). The result is 
Eq. (fl2|) in Section 12.2.21 One can analytically show that for a time t after 
photoionization 
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Ties = T i>ave + T c (0)t 2 /(r c 2 xp + t 2 ). (B.4) 
The characteristic plasma expansion time is r cxp from Eq. ff25l) . This enables 
a clear and well-defined analysis of the data, such as shown in Fig. In 
principle, it also separates the various contributions to the Doppler broadening 
of the spectrum, although the expansion energy typically becomes an order of 
magnitude greater than the thermal energy, making quantitative measurement 
of the thermal energy difficult at later times. One can also consider comparison 
of Eqs. (fT2j) and (1B.4P with data as a check of the self-similar expansion model 
(Eq. flMD). 

References jM] and [7T] numerically justify approximate expressions relating 
T i eff to T i ave , T e (0), time, and other physical parameters when the ion tem- 
perature is not assumed to be constant. 



C The Quasineutral Approximation 



The assumption of quasineutrality constitutes a central part for the theoretical 
description and the analysis of ultracold plasma experiments. The condition of 
approximately equal electron and ion densities (p e ~ p ; ) implies a well defined 
charge separation, but this fact often causes confusion and is rarely discussed 
in the literature. We provide a more detailed discussion of its physical meaning 
in this section. 

We start our discussion with the set of Vlasov equations 

df e df e e df e dip . _ , . 

at or c m e av e ar e 

at ori mi av; ari 
for the phase space distribution functions f a (r a , v a ), where m a is the electron 
(a = e) and ion (a = i) mass and ip is the total mean-field potential determined 
by the Poisson equation 

A<p , (C.3) 

with paipha = / f a dv a . Most straightforwardly, the solution of these equations 
is obtained by calculating the potential ip from Eq. (1C.3I) at every instant of 
time to propagate the phase space densities according to eqs. OC.ip and (1C2[) . 
However, the Vlasov equations could be used alternatively to calculate the 
electrostatic potential, while the Poisson equation determines the total charge 
density in this case. To do this we consider the hydrodynamic velocities 



Na 1 / v Q / a (r a ,v a )dv Q (C.4) 
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whose time evolution is obtained by multiplying eqs. (lC.lj) and (1C.2j) by v c 
and v i; respectively, and integrating over velocities. Upon subtracting of the 
resulting equations and using the continuity equations for the electron and ion 
densities we obtain for the electrostatic potential 



For isotropic, thermal velocity distributions the pressure tensors V a become 
simple scalar fields given by V a = k^T a p a . Since the time evolution of the hy- 
drodynamic velocities and the pressure tensors follows from the self-consistent 
solution of the Vlasov equations this expression is completely equivalent to 
the Poisson equation (1C.3I) . While the first term in Eq. (1C.5I) describes lo- 
cal charge imbalances between the two species, the second term arises solely 
from the thermal pressure of the particles. The latter only vanishes for the 
unlikely special case of V\ = ^r'Pe, and, hence, yields a finite contribution to 
the electrostatic potential in general. On the other hand, the first can be ne- 
glected within the quasineutral approximation, requiring almost equal charge 
densities and local currents, i.e. pi ~ p e and u ; = u c . Thus the electrostatic 
potential can be calculated from 



which coincides with Eq. (1211) . introduced in section 13.1.11 In the second step 
we have employed an adiabatic approximation for the electrons, making use 
of the relation m e -C m;. 

Note, that a direct application of the quasineutrality condition p e w p x in the 
Poisson equation (1C.3|) would have resulted in a vanishing charge separation 
potential. The preceding discussion however shows that this is not the case 
since the thermal pressures of the plasma components generally yield a finite 
contribution of zero order in the charge separation. Rather the potential is 
given by Eq. (1C.6I) . while the Poisson equation determines the separation of 
charge densities 5 = p\ — p e of the two species. However, we have 5/p e <C 1, 
ensuring quasineutrality and corresponding to the condition for the existence 
of a plasma state, introduced by Eq. (I6TI) . 
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